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The Collins-Williams Regge calculus models of FLRW space-times and Brewin’s subdivided mod¬ 
els are applied to closed vacuum A-FLRW universes. In each case, we embed the Regge Cauchy 
surfaces into 3-spheres in and consider possible measures of Cauchy surface radius that can be 
derived from the embedding. Regge equations are obtained from both global variation, where entire 
sets of identical edges get varied simultaneously, and local variation, where each edge gets varied 
individually. We explore the relationship between the two sets of solutions, the conditions under 
which the Regge Hamiltonian constraint would be a first integral of the evolution equation, the 
initial value equation for each model at its moment of time symmetry, and the performance of the 
various models. It is revealed that local variation does not generally lead to a viable Regge model. 
It is also demonstrated that the various models do satisfy their respective initial value equations. 
Finally, it is shown that the models reproduce the correct qualitative dynamics of the space-time. 
Furthermore, the approximation’s accuracy is highest when the universe is small but improves over¬ 
all as we increase the number of tetrahedra used to construct the Regge Cauchy surface. Eventually 
though, all models gradually fail to keep up with the continuum FLRW model’s expansion, with 
the models with lower numbers of tetrahedra falling away more quickly. We believe this failure to 
keep up is due to the finite resolution of the Regge Cauchy surfaces trying to approximate an ever 
expanding continuum Cauchy surface; each Regge surface has a fixed number of tetrahedra and as 
the surface being approximated gets larger, the resolution would degrade. Finally, we note that all 
Regge models end abruptly at a point when the time-like struts of the skeleton become null, though 
this end-point appears to get delayed as the number of tetrahedra is increased. 
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I. INTRODUCTION 

Since its initial formulation in 1961, Regge calculus 
[H , a discretisation of general relativity, has remained an 
enduring subject of interest, particularly because of its 
applications to numerical relativiW and to path-integral 
formulations of quantum gravity @t 3| • The focus of this 
paper is on its numerical applications. Regge calculus 
can, in principle, approximate any solution to the Ein¬ 
stein held equations of general relativity, and indeed, it 
has been shown that solutions of Regge calculus converge 
at second order in the discretisation edge-lengths to so¬ 
lutions of the Einstein held equations Q . Regge calculus 
has been applied to approximate a wide range of space- 
times [9l-ll4| , and comparison of the Regge approximation 
with the continuum solution, when known, has demon¬ 
strated good agreement, indicating that Regge calculus 
offers a viable approximation to general relativity. 

The main value of Regge calculus though is that it can 
be used to approximate space-times where exact solu¬ 
tions are difficult to obtain. In particular, it offers a non- 
perturbative approach to approximating space-times, po¬ 
tentially elucidating certain physics that, owing to the 
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non-linear nature of the Einstein held equations, would 
not be evident using perturbative approaches. Eor ex¬ 
ample, it has been argued that to unambiguously under¬ 
stand the effects of inhomogeneities on cosmological ob¬ 
servables, a non-perturbative model of inhomogeneities 
is necessary |15l - [l7| : Regge calculus may assist in this 
respect. Regge calculus can also complement other ap¬ 
proaches to numerical relativity: when no exact solution 
is available for comparison, consistency of results between 
different numerical approaches would strengthen conh- 
dence in any conclusions. 

The focus of this paper is on further exploring a Regge 
approximation to the Friedmann-Lemaitre-Robertson- 
Walker (ELRW) space-times. We shall focus on a Regge 
formalism hrst devised by Collins and Williams (CW) 
[l^ and subsequently expanded by Brewin [I^. This 
formalism was hrst applied to closed dust-hlled FLRW 
cosmologies, and the resulting approximation was able 
to reproduce the behaviour of the continuum space-time 
rather well. In this paper, we shall test the CW formal¬ 
ism against another cosmological space-time for which 
the exact solution is well-known, the closed FLRW uni¬ 
verse with non-zero cosmological constant A. A key dif¬ 
ference between this universe and that studied by Collins, 
Williams, and Brewin is that the dust-hlled universe 
eventually collapses back in on itself whereas the A- 
FLRW universe expands indehnitely. By comparing the 
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approximation against the exact solution, we can learn 
more about the CW formalism’s performance and range 
of applicability in general. As this paper will show, we 
have also uncovered several interesting properties of the 
CW formalism that were not explored in the original pa¬ 
pers. With a better understanding of this formalism, we 
should then be better-equipped to adopt it, for instance, 
to modelling inhomogeneous cosmologies where exact so¬ 
lutions are not known. As an example, drawing on the 
lessons learned here, we have recently applied this for¬ 
malism to model ‘lattice universes’ wherein the matter 
content consists of point particles distributed into a reg¬ 
ular lattice [ 2 ^. 0 

The FLRW metric is founded upon the Copernican 
principle which posits that the universe admits a foliation 
of constant-time Cauchy surfaces such that each surface 
is perfectly homogeneous and isotropic; this leads to a 
family of metrics that can be expressed in the form 


( 1 ) 

where a[t) is a time-dependent function known as the 
scale factor, t is the time parameter, and fc is a curvature 
constant. The sign of k determines whether the constant- 
t Cauchy surfaces are open, flat, or closed, with k < 0 
being open, k = 0 being flat, and fc > 0 being closed. 

By inserting this metric into the Einstein field equa¬ 
tions, one obtains a pair of differential equations for aft) 
known as the Friedmann equations; these are 
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where p and p are the energy density and pressure of any 
fluid filling the space, and A is the cosmological constant. 

For closed vacuum universes with a non-zero cosmo¬ 
logical constant, we have that k > 0, p = p = 0, and 
A ^ 0; solving the Friedmann equations then yields 


where the integration constant has been chosen so that 
d{t) = 0 when t = 0, and where we have chosen a scaling 
of aft) such that it corresponds to the radius of curva¬ 
ture of constant-f Cauchy surfaces; for such a scaling, the 
curvature constant fc would in turn be re-scaled to fc = 1. 


^ There has been great interest recently in studying such toy uni¬ 
verses to try and understand whether dark energy is indeed 
needed to explain the observed acceleration of the universe or 
whether this acceleration can be explained away as an apparent 
effect arising from the influence of the late universe’s inhomoge¬ 
neous matter distribution on cosmological observables [ 21 II . 


Constant-f Cauchy surfaces of any closed FLRW 
universe can always be embedded as 3-spheres in 4- 
dimensional Euclidean space Such an embedding 

requires scaling fc to be unity so that aft) equals the 
3-sphere radius. The embedding is then given by 

r = siny, 

= aft) cosy, 

= a(t) sinycosS, (5) 

= aft) sin y sin 9 cos </), 
x^ = aft) sin x sin 9 sin (f, 

for 0 < X, 0 < TT and 0 < (j) < 27r. Such a 3-sphere would 
have a volume of 


C^FLRw(i) = 2Tr^a(t)^, (6) 

and an expansion rate of 

t^FLRw(i) = 67r^a(f)^ d{t). (7) 

The FLRW metric can then be expressed as 

= —dt^ + a‘^{t) [dy^ -|- sin^ y (dd^ -|- sin^ dd(/)^)] . 

( 8 ) 

Regge calculus [l| approximates any curved space-time 
using a piece-wise linear manifold constructed out of flat 
4-blocks: the blocks are ‘glued’ together such that neigh¬ 
bouring blocks share an entire 3-face, and as the blocks 
are flat, the metric inside is the Minkowski metric. Regge 
space-times are generally referred to as skeletons. Curva¬ 
ture in a skeleton manifests itself as conical singularities 
concentrated on the sub-faces of co-dimension 2; these 
are known as hinges. If a hinge were flat, then the dihe¬ 
dral angles between all 3-faces meeting at it would sum 
to 27r; any deviation from 27r provides a measure of the 
curvature and is known as the deficit angle. The edge- 
lengths serve as the Regge analogue of the metric and are 
determined by the Regge field equations, a set of equa¬ 
tions analogous to the Einstein field equations of general 
relativity; and just as the Einstein field equations can be 
obtained by varying the Einstein-Hilbert action with re¬ 
spect to the metric, so can the Regge field equations be 
obtained by varying a Regge action with respect to the 
edges. 

One of the central features of the CW formalism is its 
skeleton, designed to strongly mirror the structure and 
symmetries of the continuum FLRW space-time it is ap¬ 
proximating M- In analogy with FLRW space-times, 
CW skeletons are also foliated by a one-parameter fam¬ 
ily of space-like Cauchy surfaces. But now, each surface 
is essentially a triangulation of an FLRW 3-sphere using 
equilateral tetrahedra such that all vertices, edges, and 
triangles are identical to each other; in this way, the sur¬ 
faces would mimic as closely as possible the Copernican 
symmetries. According to Coxeter [ 2 ^, such a triangula¬ 
tion of the 3-sphere is only possible with 5, 16, and 600 
tetrahedra; Table |T] summarises the numbers of vertices. 
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Table 1. The number of simplices in each of the three triangu¬ 
lations of the 3-sphere with equilateral tetrahedra as well as 
the number of triangles meeting at any edge. We introduce 
Ns, N 2 , Ni, and Nq to denote the numbers of parent tetra¬ 
hedra, triangles, edges, and vertices in the Cauchy surface. 


Tetrahedra 

m 

Triangles 

(A2) 

Edges 

(fVi) 

Vertices 

(No) 

Triangles 
per edge 

5 

10 

10 

5 

3 

16 

32 

24 

8 

4 

600 

1200 

720 

120 

5 


edges, triangles, and tetrahedra for each case. As with 
FLRW Cauchy surfaces, all CW Cauchy surfaces are re¬ 
quired to be identical to each other apart from an overall 
scaling, represented by the length l{ti) of the tetrahedral 
edge, ti being a discrete time parameter labelling the foli¬ 
ation; thus l{ti) would be a Regge analogue of the FLRW 
scale factor a{t). To complete the skeleton’s construction, 
the CW Cauchy surfaces are glued together by a series of 
time-like edges connecting vertices in one surface to their 
time-evolved images in the next; these edges are known 
as struts. The world-tubes of the tetrahedra between two 
consecutive Cauchy surfaces would then correspond to 
the 4-blocks for this skeleton. To ensure homogeneity of 
all vertices in any Cauchy surface, all struts between pairs 
of consecutive Cauchy surfaces are required to be identi¬ 
cal to each other. Finally, by taking the limit where the 
separation between surfaces goes to zero, one can gen¬ 
erate a continuum time formulation of Regge calculus. 
Collins and Williams first applied their construction to 
model closed dust-filled FLRW universes, and the con¬ 
tinuum time function l{t) for the edge-lengths behaved 
very similarly to the equivalent FLRW scale-factor a{t), 
with models with a greater number of tetrahedra yielding 
better accuracy. 

The CW formulation was further explored and ex¬ 
tended by Brewin M- In particular, Brewin devised an 
algorithm to triangulate each tetrahedron into smaller 
tetrahedra thereby generating secondary models. The 
virtue of this algorithm is that it can in principle be re¬ 
peated indefinitely, thereby yielding even finer approx¬ 
imations to the underlying FLRW surfaces. However, 
subdivision comes at the expense of some of the sym¬ 
metries inherent in the original CW surfaces: the new 
tetrahedra would no longer be identical nor necessarily 
equilateral. For instance, after the first generation of sub¬ 
division, each Cauchy surface would instead have three 
sets of vertices, three sets of edge-lengths, three sets of 
triangles, and three sets of tetrahedra. All members of 
any set would be identical to each other, and in this sense, 
there was still some Copernican symmetry. We shall re¬ 
fer to CW’s original models as the parent models and any 
subdivided ones as children models. The children mod¬ 
els were also applied to closed dust-filled FLRW universes 
and were found to better approximate the continuum uni¬ 
verse compared to their parent counterparts, again with 


accuracy commensurate with the number of tetrahedra. 

Both parent and children Cauchy surfaces can also 
be embedded in 3-spheres in E^, effectively embedding 
them in FLRW Cauchy surfaces like those described by 
dSD- Such an embedding was first mentioned by Collins 
and Williams and more extensively developed by Brewin. 
Since the CW Cauchy surfaces are approximations of 
3-spheres, the embedding would map all vertices from 
a parent Cauchy surface into a single 3-sphere, and it 
would, in general, map each of the three sets of vertices 
in a child Cauchy surface into its own 3-sphere. However 
in Brewin’s embedding of the child Cauchy surface, he 
considered only the case where the three sets of vertices 
had been constrained to lie on the same 3-sphere. He then 
proposed that, for both the parent and children models, 
the corresponding 3-sphere radius R{t) could provide a 
more accurate analogue to the FLRW scale factor a{t) 
than the tetrahedral edge-length. In this paper, we shall 
generalise Brewin’s investigation in two ways: we shall 
look at alternative radii that can also serve as analogues 
to a{t), and we shall also consider models where the three 
sets of children vertices do not necessarily lie on the same 
3-sphere. 

Brewin has also pointed out certain analogies between 
the ADM formalism and the CW formalism. He has 
likened the tetrahedral edge-lengths to the 3-metric of 
an ADM foliation and the Regge equations obtained from 
varying the tetrahedral edges to the evolution equations. 
He has also likened the struts and diagonals to the ADM 
lapse and shift functions, respectively, and the Regge 
equations obtained from their variation to the Hamilto¬ 
nian and momentum constraints, respectively. Thus in 
this paper, we shall refer to the Regge equations obtained 
from the tetrahedral edges as the evolution equations, 
from the struts as the Hamiltonian constraints, and from 
the diagonals as the momentum constraints. 

This paper is organised as follows. The second section 
provides a brief review of Regge calculus in general and 
the CW formalism specifically. The third section applies 
the parent CW formalism to the A-FLRW universes. We 
begin with an exposition of the 4-block’s geometry. We 
follow this by presenting possible embeddings of parent 
CW Cauchy surfaces into 3-spheres and discussing pos¬ 
sible measures for the Cauchy surfaces’ radii that can 
serve as analogues to a{t). We then vary the Regge ac¬ 
tion in two different manners to arrive at the Regge held 
equations. In the hrst approach, we impose the sym¬ 
metry constraints on the edges hrst before varying the 
action. The skeleton must then continue satisfying the 
constraints even under variation; so if we vary one edge, 
then all edges constrained to share the same length must 
get varied at the same time: we would simultaneously 
vary either all tetrahedral edges in a Cauchy surface or 
all struts between a pair of consecutive surfaces. This 
was the approach followed by Collins and Williams, and 
we shall call it global variation. From this approach, we 
hnd that the CW Hamiltonian constraint for our model 
is actually a hrst integral of the evolution equation; thus. 
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the Hamiltonian constraint is sufficient to determine our 
model’s evolution. In the second approach, we instead 
vary each edge individually first and afterwards impose 
the symmetry constraints on the resulting Regge equa¬ 
tions; this is the more standard way of varying the action 
in Regge calculus, and we shall call this local variation. 
In this case, we find that the Hamiltonian constraints 
are again first integrals of the evolution equations pro¬ 
vided we also satisfy the momentum constraints. How¬ 
ever, these momentum constraints impose rather unphys¬ 
ical constraints, and we therefore dismiss the local model 
as unviable. The global and local Regge equations can 
also be related through a chain rule; when we do this, 
we again arrive at the same conclusion about the local 
model’s unviability. We next consider the initial value 
equation in the context of Regge calculus and demon¬ 
strate that the global models do satisfy this equation at 
its moment of time symmetry. We conclude our investi¬ 
gation of the parent models with a brief speculation on 
the reasons for the local models’ breakdown before finally 
examining the evolution of the global models. 

In the final section, we turn our attention to the sec¬ 
ondary models obtained by subdividing the tetrahedra 
of the parent models. We begin by presenting Brewin’s 
scheme for subdividing the CW models. This is followed 
by the embedding of children Cauchy surfaces into 3- 
spheres. We next present the geometric quantities needed 
to compute the varied Regge action. After briefly dis¬ 
cussing the local variation of the action, we determine 
the global Regge equations obtained by varying with re¬ 
spect to both the struts and the tetrahedral edge-lengths, 
and we briefly discuss the conditions under which the 
Hamiltonian constraint would be a first integral of the 
evolution equation. We then consider the initial value 
equation for the children models and demonstrate that 
the models do satisfy this equation. Whereas our method 
for applying the equation to the parent models is based 
on a proposal by Wheeler (2^, we have devised a differ¬ 
ent form of the equation for the children models; this is 
because Wheeler’s proposal can only be applied to Regge 
Cauchy surfaces where there is a well-defined ‘volume per 
vertex’, which is the case for the parent models but not 
for the children models. To the best of our knowledge, 
our alternative form of the initial value equation is com¬ 
pletely novel and can be applied to any time-symmetric 
Regge Cauchy surface in general. Finally, we examine 
the evolution of the children models, comparing their 
performance against those of the parent models, and we 
speculate briefly on how the children models might be 
extended. 

In this paper, we shall use geometric units where G = 
c = 1 . 


II. REGGE GALGULUS AND THE CW 
SKELETONS 

In general relativity, the Einstein field equations can 
be derived by varying the Einstein-Hilbert action 

Seh = J (R - 2 A) d^x, (9) 

with respect to the metric tensor 5 ^ 1 ,, where A is the 
cosmological constant, R is the Ricci scalar, and g = 
det{g^,y). 

When applied to a piece-wise linear manifold, as found 
in Regge calculus, this reduces to the Regge action [ij 

SRe,,e = ^ I E - E j ’ (10) 

\i G {hinges} i G {4-blocks} j 

where Ai is the area of a hinge in the Regge skeleton. Si 
its corresponding deficit angle, and 14^^^ the volume of 
a 4-block. The first summation is over all hinges in the 
skeleton while the second is over all 4-blocks. The deficit 
angle Si at hinge i is given by 

5, = 27r-E^f, (11) 

J 

where 0^0 is the dihedral angle between the two faces of 
block j meeting at the hinge and the summation is over 
all blocks meeting at the hinge. 

Since the skeletal edge-lengths are the Regge analogue 
of the metric, the Regge action is varied with respect to 
an edge-length £j to get the Regge field equation 



where the variation of the deficit angles has cancelled out 
owing to the well-known Schlaffi identity [I|, 


this identity holds for any individual block fc, with the 

(i) 

summation being over all hinges in the block and 0^ 
being the block’s dihedral angle at hinge Ai^ 

Regge calculus customarily uses simplicial manifolds 
where the skeleton’s fundamental building block is the 



^ In the standard formulation of Regge calculus, one actually uses 
a simplicial manifold where every block is a 4-simplex, and the 
Schlafli identity is usually formulated in terms of simplices rather 
than arbitrary blocks. However, any block can always be trian¬ 
gulated into simplices, and one can then apply the simplicial 
form of the Schlafli identity to the triangulated block to obtain 
the form of the identity we have above, using the chain rule if 
necessary to satisfy any constraints on the block’s geometry. 
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4-simplex. The geometry of a single n-simplex can al¬ 
ways be completely determined by specifying the lengths 
of its C(n -I- 1, 2) edges; therefore, the geometry of an en¬ 
tire simplicial skeleton can be completely determined by 
specifying the lengths of all its edges. However, the fun¬ 
damental building blocks of CW skeletons are instead 4- 
blocks corresponding to the truncated world-tubes of the 
tetrahedra as they evolve from one Cauchy surface to the 
next, as illustrated in Fig. [TJ Without extra constraints 
on the 4-block’s internal geometry, the skeleton’s geom¬ 
etry would not be completely determined. As an anal¬ 
ogous 2-dimensional example, consider the geometry of 
a quadrilateral where only its four external edge-lengths 
are known: there is a wide range of possible quadrilat¬ 
erals that would have these four edge-lengths, ranging 
from trapezia to irregular quadrilaterals; without extra 
constraints, it would be impossible to determine a unique 
quadrilateral. 


D' 



Figure 1. An equilateral tetrahedron of edge-length h at time 
ti evolves to a tetrahedron of edge-length h+i at time ti+i, 
tracing out a 4-dimensional world-tube. The struts are all of 
equal length. 

The CW 4-block’s geometry is therefore constrained by 
requiring the tetrahedron to simply expand or contract 
uniformly about its centre as it evolves from one end 
of the 4-block to the other. Brewin has indicated that 
such a geometry is equivalent to imposing the following 
two requirements on the 4-block: (i) that all struts have 
the same length; and (ii) that there be no twist or shear 
along the 4-block. He has likened these requirements to 
a choice of lapse and shift function in the ADM formal¬ 
ism. Indeed, the standard form of the FLRW metric © 
also implies a certain foliation of FLRW space-time, and 
Collins and Williams’ choice seems closest to the lapse 
and shift implicit in this foliation. 

The 4-block geometry can be described by introducing 
a Cartesian co-ordinate system into the 4-block. We shall 
denote by the Cauchy surface at time U. Let the 
tetrahedron in have length li := and label its 


vertices by A, B, C, and D. The origin is then set to 
be at the tetrahedron’s centre; the a;-axis is taken to lie 
parallel to AB, the y-axis to pass through vertex C, and 
the z-axis to pass through vertex D. The vertices’ co¬ 
ordinates are therefore^ 


A 

B 

C 

D 


li 1% 


d. 


h 


2’ 2^3’ 2^/6' 
li li li 

2’ 2VE 

„ li 

’ 73 ’ 276 ’^^ 


272 ’ 


(13) 


To help simplify future calculations with this geome¬ 
try, we have chosen to use a Euclidean metric; we have 
therefore introduced the imaginary unit l in our time co¬ 
ordinate so that inner products would effectively yield a 
signature of (+, +, +, -). 

This tetrahedron will evolve to another of edge-length 
lij^i := l{tij^i) in surface We shall denote the im¬ 

age of vertices A, B, C, D in by A', B', C, D', 
respectively. Their co-ordinates are given by an analo¬ 
gous expression to uni), where each vertex is replaced by 
its primed counterpart and each subscript i by i-1-1. This 
ensures the struts all have the same length. Additionally, 
in this co-ordinate system, the tetrahedron would simply 
expand or contract uniformly about its centre in the spa¬ 
tial dimensions, as required. Figure [T] illustrates this par¬ 
ticular 4-block. For simplicity, we shall sometimes refer 
to the tetrahedron in as the upper tetrahedron and 
the one in as the lower tetrahedron, as that is how 
they appear in the figure. 

The CW constraints on the 4-block geometry are im¬ 
posed in different manners depending on whether we are 
globally or locally varying the skeleton. Under local vari¬ 
ation, we must first completely triangulate the skeleton, 
thereby generating a completely simplicial manifold. All 
edges in this fully triangulated skeleton are then consid¬ 
ered independent of all others, including the newly intro¬ 
duced edges; when one edge is locally varied, all others 
are held constant. After the Regge equations are ob¬ 
tained, the constraints on the geometry are then imposed; 
this is done by setting the various edges to have the ap¬ 
propriate lengths consistent with the CW 4-block geom¬ 
etry: that is, all tetrahedral edge-lengths in a Cauchy 
surface would be set equal; all strut-lengths between a 
pair of consecutive surfaces would be set equal; and all 
diagonal-lengths between a pair of consecutive surfaces 
would be set equal. 


® We could more generally have used a time co-ordinate iTi := 
LT{ti) instead for the co-ordinates in II13II . but this would still 
lead to the same set of final equations. So for simplicity, we 
choose to use ti as the vertices’ time co-ordinates. 
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The specific 4-block represented by Fig. [T] and co¬ 
ordinate system (fT51) is triangulated by introducing di¬ 
agonals AD', BD', CD', AC, BC, and AB'-, this cor¬ 
responds to one diagonal above each of the tetrahedral 
edges, as illustrated in Fig. [2l and these diagonals di¬ 
vide the 4-block into four distinct 4-simplices, ABCDD', 
ABCC'D', ABB'C'D', and AA'B'C'D'. One can gen¬ 
erate a consistent triangulation of the entire skeleton by 
triangulating each of its 4-blocks in this manner. The 



Figure 2. The world-sheets generated by the six tetrahedral 
edges and their triangulation into triangular time-like hinges. 

CW geometry is then imposed on a 4-block by requiring 
its lower tetrahedral edges to have length k, its upper 
tetrahedral edges to have length k+i, its struts to have 
length nii given by 

"1? = - l) (14) 


where we have introduced the notation 



and its diagonals to have length di given by 

di = ^ {^k+i + li)^ — . (15) 

These lengths are all derived from the co-ordinates m 
and its counterpart. 

Under global variation, the 4-block geometry is instead 
constrained to have the CW geometry before the skele¬ 
ton gets varied, and this is done without introducing any 
new edges. Global variation preserves the CW geometry 
because when a tetrahedral edge-length gets varied, all 
tetrahedral edges in the same Cauchy surface gets var¬ 
ied; in each 4-block, the varied tetrahedron simply gets 
re-scaled uniformly about its centre. Similarly, when a 
strut-length gets varied, all struts between the same pair 
of surfaces gets varied. In either case, each 4-block sim¬ 
ply becomes a different CW 4-block. We note that under 
global variation, each Cauchy surface would consist of 


only two independent edges, the tetrahedral edge and 
the strut. 

As we mentioned in the Introduction, local variation 
is the standard way of doing Regge calculus; it is more 
similar to how standard general relativity is done. In 
general relativity, global variation would be analogous to 
requiring the metric in the Einstein-Hilbert action to be 
of FLRW form ([T|) , and then varying the action with re¬ 
spect to a(t); this imposes the Copernican symmetries 
prior to varying the Einstein-Hilbert action. The stan¬ 
dard approach would be to vary the Einstein-Hilbert ac¬ 
tion first, yielding the Einstein field equations, and then 
setting the metric to be of ELRW form. 

B rewin 0 has explored to greater depth the relation¬ 
ship between solutions of the global Regge equations and 
solutions of the local Regge equations. He has shown 
that, in general, global and local equations would not 
necessarily lead to the same set of solutions; rather the 
local solutions would, under certain circumstances, form 
a subset of the global solutions. We shall demonstrate 
this relationship explicitly for the parent model later on. 


III. PARENT MODELS OF CLOSED 
VACUUM A-FLRW UNIVERSES 

A. Embedding Cauchy surfaces into a 3-sphere 

As Collins and Williams first noted and Brewin fully 
explored, a CW Cauchy surface can be embedded into 
E^ such that all vertices lie on a 3-sphere of radius Ri. 
It is most natural to parametrise this 3-sphere using a 
set of polar co-ordinates (x,0,^). If (x^, , x'^) is a 

set of Cartesian co-ordinates in E^, then points on the 
3-sphere are can be parametrised by 

= Ri cosy, 
x^ = Ri sin Y cos 0, 

, (16) 
X = Ri sin X sin 9 cos (j), 

= Ri sin X sin ^ sin (j), 

with X, 0 £ [0,7r] and </> G [0, 27r). This embedding pro¬ 
vides a natural framework within which to study and 
elucidate the underlying geometry of the CW Cauchy 
surface. Most importantly, it makes clearer the relation¬ 
ship between the CW Cauchy surfaces and the FLRW 
3-spheres they approximate. 

We begin by embedding a subset of the vertices into 
the 3-sphere; a schematic diagram of our embedding is 
given in Fig. [31 By symmetry of the 3-sphere, we can 
always choose polar co-ordinates such that one vertex is 
located at (x, ^*, </>) = (0,0,0). Then, we also have free¬ 
dom to choose 9 and (j> co-ordinates such that one of the 
neighbouring vertices is at {x,9,(j)) = (xo,0,0) for some 
Xo- If is the number of triangles meeting at an edge, 
then there will be n vertices surrounding the edge formed 
by the first two vertices. To see this, we can consider the 
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Figure 3. A schematic diagram of an equilateral tetrahedron 
of edge-length k embedded into a 3-sphere of radius Ri. One 
dimension has been projected out. 


Table II. The polar co-ordinates of n+2 neighbouring vertices. 


Vertex 

X 

e 

(p 

1 

0 

0 

0 

2 

Xo 

0 

0 

3 

XI 

00 

0 

4 

XI 

00 

27r 

n 

n + 2 

XI 

00 

2(n—l)7r 
n 


first two vertices as forming a common base for the n 
triangles and each of the n vertices as forming the apex 
of each of the n triangles. The values for n correspond¬ 
ing to the different models are listed in the final column 
of Table H We again have freedom to choose the (p co¬ 
ordinate such that one of these n vertices is located at 
(x, 9, (p) = (xi: ^* 0 ,0) for some xi and Oq. Once this choice 
is made, the remaining n — 1 vertices would be located at 
(Xi,^'o,27r/n), (xi, 6 * 0 ,47r/n), ..., (xi, 6 * 0 ,2(n - l)7r/n). 
The co-ordinates of these vertices have been summarised 
in Table im By requiring all distances to be k between any 
pair of neighbouring vertices, we obtain the equations 


We then solve these to obtain 


Xo = Xi, 
cos Xo = 


cos ^ 


COSf/o = 


1 - 2 cos ^ 

n 

cos^ 

1 - cos ^ ’ 


Zo := ^ = 72 
■th'i 


1 - 3cos^ 
_ 22^ 

1 - 2 cos ^ 


1/2 


(18) 

(19) 

( 20 ) 
( 21 ) 


These are the relations Brewin obtained for the par¬ 
ent models’ embedding, although we have derived them 
here independently of him. 

From this embedding, we see that as the edge-lengths 
l{ti) expand and contract, the 3-sphere simply expands 
and contracts about its centre, and the vertices simply 
move radially inwards or outwards accordingly. The an¬ 
gular positions of these vertices remain constant. 

Two interesting features of the Cauchy surface geom¬ 
etry come to light from this embedding. First, we can 
now see that li and Ri are related by the constant ra¬ 
tio Zq. Most notably, this ratio is independent of the 
label i and hence of time ti. Thus, we can define a ra¬ 
dius R{ti) = l{ti)/ZQ for our CW Cauchy surfaces, and 
this serves as a natural analogue to the FLRW scale fac¬ 
tor a{t). Secondly, we see that as the edge-lengths l{ti) 
expand and contract, the 3-sphere simply expands and 
contracts about its centre, and the vertices simply move 
radially inwards or outwards accordingly; their angular 
positions remain constant. 

Our embedding above has yielded one possible defini¬ 
tion of radius for the CW Cauchy surface, namely the ver¬ 
tices’ embedding radius. This was the definition Brewin 
chose as his Regge analogue to the FLRW scale factor 
of a{t). However, there are other equally plausible def¬ 
initions of radius for the Cauchy surface: we could just 
as well have chosen the radius of any other point in the 
tetrahedra, as there will always be a set of points in the 
Cauchy surface sharing that same radius. Some possibil¬ 
ities are the radius Ri to the centres of edges 


Ri{ti) 


72 U-2cosf 


1/2 

Riti), 


Ri{U) V 1 - cos ^ ) 


( 22 ) 

(23) 


(ii 2 )'=if = 2i?7l-cosxo) 

= 2-R71 - cosxi) 

{L 2 iif = Ij = 2Rj{l - cosxocosxi - sinxosinxi cos 6 >o) 


{LmY = sin^ xi sin^ 9q [l — cos — ^ , 

(17) 

where Lij denotes the distance between vertices i and j. 


the radius i ?2 to the centres of triangles 


i?2(ti) 


1{U) 

R2{ti) 


(- 


73 VI - 2 


cos' 


1/2 


R{ti), 


76(1-3 


cos 


27r\ 
n ) 


1/2 


(24) 

(25) 




















or the radius i ?3 to the centres of tetrahedra 


Rsiti) 


1{U) 

RsiU) 


1 / 1 + cos ^ 

2 l^l-2cos^ 


1/2 

R(ti), 


2 V 2 


/l-3cos^ 


1/2 


(26) 

(27) 


The main point to notice is that regardless of which ra¬ 
dius we choose, the ratio between that radius and l{ti) is 
always a constant independent of ti. 

Since no particular choice of radius seems preferred, 
one natural choice would be to average over all radii 
across the entire Cauchy surface. We have numerically 
computed this average radius R in terms of R{ti) to be 


R{ti) = < 


0.484066 i?(t,) 
0.627392 i?(t,) 
0.940901 i?(t,) 


5 tetrahedra model, 

16 tetrahedra model, 
600 tetrahedra model. 


(28) 

The derivation of these numbers has been explained in 
Appendix O 

We can also consider the effective radius R{ti) obtained 
by treating the volume of the Cauchy surface as if it were 
the volume of a 3-sphere. The volume of a 3-sphere of 
radius R{t) is 


Uit) = 2n‘^R{tf, (29) 

while the volume of an A^s-tetrahedra CW universe is 

Nn 

UNAU) = -^l{Uf. (30) 

If we equate the two expressions, we find that the effective 
radius is 


Finally, we shall consider one more possible definition 
of the radius when we have obtained the equations for 
l{t) and l{t) in the continuum time limit. Like all other 
radii, this radius R{t) is related to l{t) by a constant Z, 

R{t) = Zl{t), (34) 

and we define Z by requiring R = a{t) when both 
dR/dt = 0 and a = 0. 


B. Global variation of the parent models 

Under the global constraints, the CW skeleton simpli¬ 
fies significantly. Most importantly, there are now only 
two distinct sets of hinges associated with any Cauchy 
surface. The first corresponds to the world-sheets of the 
tetrahedral edges between surfaces and E^+i. These 
hinges are time-like and trapezoidal; an example would 
be hinge ABA'B' in Fig. [T] The second corresponds to 
the equilateral triangular faces of the tetrahedra. These 
hinges are space-like and have edges of length k; an ex¬ 
ample would be the hinge ABC in Fig. [TJ 

Thus for a skeleton satisfying the global constraints, 
the Regge action can be written as 

4-blocks} 

(35) 

We shall now derive the Regge equations by global vari¬ 
ation. As mentioned previously, there are only two dis¬ 
tinct types of edges characterising a global skeleton, the 
struts and the tetrahedral edges. If we vary the action 
with respect to a strut, we obtain the Hamiltonian con¬ 
straint 


/ /V 3 \ 

For comparison with the average radii in (1281) , the nu¬ 
merical values for R{t) in terms of the vertex radius R{ti) 
are 


R{ti) = < 


0.490488 R(t0 
0.646482 R{U) 
0.945651 R{U) 


5 tetrahedra model, 

16 tetrahedra model, 

600 tetrahedra model, 

(32) 


and the fractional difference between these numerical fac¬ 
tors and those in (E51) are 


Ritj) - R(ti) 
R{ti) 


-0.0132668 

-0.0304271 

-0.00504837 


5 tetrahedra model, 

16 tetrahedra model, 
600 tetrahedra model. 


j33) 

The two radii are very close to each other, but with R{ti) 
consistently greater than R{ti) by a slight amount. 




drui 


drrii 


(36) 


where Ni and A 3 are the numbers of edges and tetrahe¬ 
dra in a Cauchy surface and equal the numbers of trape¬ 
zoidal hinges and 4-blocks, respectively, between any two 
consecutive Cauchy surfaces E^ and E^+i. If we vary with 
respect to a tetrahedral edge, we obtain the evolution 
equation 


Ni 


I rtrap 


dk 


= N 3 A 


dk 

dV, 


f) Atri 

AT ;:tri 


(4) 


(37) 


dk dk 


where N 2 is the number of triangles in a Cauchy sur¬ 
face; all three numbers Ni, N 2 , and N 3 are given in Ta¬ 
ble H From these equations, we see that there are only 
three types of geometric quantities relevant to the global 
Regge equations: the varied hinge areas, the correspond¬ 
ing deficit angles, and the varied 4-volumes. We shall 
now derive each in turn. 
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The area of any trapezoidal hinge between and E^+i 
is 




-I 1/2 


-^ih+i - h)'^ - ml 


(38) 


while the area of any triangular hinge in is 





(39) 


If the two hinge areas are varied with respect to mj , only 
the variation of will be non-zero, yielding 


drrij 


~ + h) 



kf-'^ 


- 1 - 1/2 


(40) 

If the space-like triangular hinges are varied with respect 
to li, we obtain 


' _ V3 
dk ~ ~ 


(41) 


If the trapezoidal hinges are varied, there will actually be 
two sets of hinges that get affected because each edge k is 
attached to two trapezoidal hinges, one between surfaces 
Ei and E^+i and the other between E^ and Ei_i. Varying 
a ‘future’ hinge with respect to k yields 


g^trap 

dh 


I' + l 4” ^2 



(42) 


and varying the corresponding ‘past’ hinge yields 


_ i \kili - k-i) - 
2 


(43) 


In general, the deficit angle for any hinge would be 
given by dm- But because all simplices are identical, 
the deficit angle on a trapezoidal hinge can be simplified 
to 


= 27r-n6»i, (44) 


where n is the number of faces meeting at the hinge and 
9i is the dihedral angle between any two adjacent faces. 
Since each trapezoidal hinge corresponds to the world- 
sheet of a tetrahedral edge and each face on this hinge 
to the world-tube of a triangle at this edge, n is equal to 
the number of triangles meeting at an edge; this number 
is listed in the last column of Table HI 

Faces ABCA'B'C' and ABDA'B'D' meeting at hinge 
ABA'B' will be separated by a dihedral angle of 0^; 
hence, 9i can be determined from the scalar product 
of the two faces’ unit normals. Let rii denote the unit 
normal pointing into ABCA'B'C and n 2 the unit nor¬ 
mal out of ABDA'B'D'] then in co-ordinate system dT51) . 
they have components 


— 


(o. 




(i-i'/) 


1/2 


and 


rin = 


(0,-2V2,1v^4) 

/ . \ 1/2 ’ 
3 ( 1-21 


and therefore 9i is given by 


cos 9i 


1 + 14^ 


(46) 


(47) 


Four faces will meet at a triangular hinge in E^. For 
hinge ABC in Fig. [H three of these faces are ABCD, 
ABCA^B'^C\ and ABCA^B^C^, where we use super¬ 
scripts 1 and t to denote the counterparts to vertices A, 
B, C in Ei_|_i and Ei_i, respectively. The fourth face 
corresponds to the neighbouring tetrahedron, which we 
denote by ABCE. By symmetry, ABCD and ABCE 
will form the same dihedral angle with ABCA^ B^ and 
with ABCA^B^C^; hence, there will only be two distinct 
dihedral angles surrounding this hinge. We take (j)\ to be 
the angle ABCD and ABCE form with ABCA^ B^ , 
and (j)\ to be the angle they form with ABCA^B^C^. 
Thus the deficit angle of ABC is 

5“ = 27r-2<^1-2<(.l. (48) 


We just used vectors orthogonal to faces to calculate 
the previous dihedral angle, but with this approach, there 
may be uncertainties over the vectors’ correct relative 
sign. So we shall take a slightly different approach here; 
we shall instead use unit vectors tangent to the two faces 
but orthogonal to the hinge. Since the hinge has co¬ 
dimension 2, there will always be a unique tangent vector 
satisfying these constraints for each face. Our approach 
is depicted schematically in Fig. SI 

We shall deduce (j)\ first. In the 4-block co-ordinates of 
m, the unit vector tangent to ABCD but orthogonal 
to ABC is simply 


'^ABCD — (Oj Oi 1) 0)j 

while the equivalent vector tangent to ABCA^ B^ is 

'“aBCAIBICI — 

1 

Thus (pl is given by 

cos((>| =-(49) 

sJlJ! -1 

By swapping for U-i, which appears implicitly in 
li^ we immediately obtain the corresponding expression 
for (j)\, that is, 



cos (pi = 


1 i 
2 V 6 

24b-l ^ 


(45) 


(50) 
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Figure 4. Two faces, separated by a dihedral angle of 
meet at a hinge. The system has been projected onto the 
plane orthogonal to the hinge. The deficit angle can be com¬ 
puted by taking the scalar product of the two vectors orthog¬ 
onal to the faces, ni and n 2 , or by taking the product of the 
two vectors tangent to the faces but orthogonal to the hinge, 
ui and U 2 . However, if we use the orthogonal vectors, we 
must be careful about their relative sign, otherwise we may 
end up taking the product with the incorrect vector, as ex¬ 
emplified by —n 2 . Instead, we have no such uncertainty if we 
work with the tangent vectors. 


The final geometric quantities required for the Regge 
equations are the 4-block volumes. The volume of 
ABCDA'B'C'D' is given by 


yA) = t;2 , ;2 




24^2 

Varying this with respect to rrij yields 




1/2 


(51) 


tween Ej and Ej+i. Varying a ‘future’ 4-block yields 


dl, '^*'24^2 


■ {h+l — h)+’<TT'i 


and varying a ‘past’ 4-block yields 


(53) 


dl, *'24^2 



— (^f i+2/i i^i + 3 Z^) 

(54) 


We can now substitute these geometric quantities into 
the Regge equations above. For the moment, we shall 
only do this for the Hamiltonian constraint (1361) . which 
yields 


= 12^2 


Ni 

N 3 A 


\j{h+i - h)"^ - 


1/2 

(27r 


n9i). 


. 

The other equation simplifies greatly in the continuum 
time limit, so we shall only present its continuum time 
form later on. 

We now take the continuum time limit, where 6 ti —>■ 0, 
to obtain a differential equation for l{t). In this limit, 
the tetrahedral edge-lengths and dihedral angles take the 
form 


dV, 


(4) 


dm, ~ 24 ^ 2 ""*^^*+^^^*^^^*+'^^*^ 


-Xk+I - liY - m] 


-1-1/2 


5ij. 


(52) 


When entire 4-blocks are varied with respect to the tetra¬ 
hedral edge-lengths, the situation is similar to the trape¬ 
zoidal hinges: each edge Ij is associated with a ‘past’ 
4-block between Ej and Ej_i and a ‘future’ 4-block be- 


h —^ 

I dt dt^ 0{diX^ , 

li — \ —^ I (t) — i dt ~b 2 ^ Oi^dt^^ , 
e, ^ 01 (t) -b0l dt + 0 {dt^) , 

0 *_i ^ eXt) + e^dt + o{dt^) , 

4)\ (X (t) + dt + 0[dt^) , 

(j)\ (^l(t) + dt + 0{dt^) , 


where in the context of continuum time, the overdot de¬ 
notes a time-derivative. The various geometric quantities 
now become 


J 


ii —^ i dt ~b ^ ^ dttX -\- o{dt ^^, 
ml — nX dt^ + 0(^dt^) , 

2 f ^ '2 

where m := { - I — 1 
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dntii 

^^trap 

dk 

dk 


I rh [ — — 1 


L 

~2 

b 

2 ~ 


- 1/2 


(i -1'") 


1/2 


+ 0 {dt ), 


1 • dt 

-II H-r- , 

2 i-izn4 


^/[+ ip _ A/4_ ^ 
2 64 


(i - iP) 


1/2 


dAf , 73, 
dk 2 ’ 


COS0, —)► 


cos 0i_i 


i+ii! 

3-ip 

i+ii! 

3 - ip 


+ 


^ di + 0(c?t^), 


(3-1^'^)' 


■ dt + 0 [dt^) , 


/l + i/2\ 

0^ = 0 ^ = arccos I -f-:- I + 0 (dt). 

V-¥V 


=- 0 ^ = - 


cos 


ll 


iV2f, + 


(i - 


2 V 6 


(i - tA) 


1/2 


^ ■ ^Idt 
I + 2 


aI 

COS®; —— 


2v^ (l - ip) 


1/2 


1_LP 

^ 24‘ , 


'i_ lidt 


= TT — (j)^ = 


arccos 


2\/6 


\(^ 


1/2 


1 - i/2 

24 . 


+ 0((ii), 


+ o(dt"), 




5K 


(4) 


drrii 


/M ? /■' - 1 


6^2 V8 


1/2 


+ 0 (dt), 


dV 


(4) 


dk 

dk 


24^2 I 2 
24^2 I 


+ dt 


-M+QI^ [-1^-1 


3.-, 


1 • dt 1 

a 

i/2 - 



-ll + 

2 1 - |P ' 

^4 

2 

64 J 


'' ^-fi-dt 
2 


+ ( -P -1 

4 


- 0(dP) 
+ 0(dP) . 


+ 0 {dP) , 
0 {dP) , 


I- 

Since 6 ^ = 9^, we shall henceforth denote this quantity tions (1551) and (1571) . we obtain, to leading order in dt, 
simply by 9. We can invert the continuum time expres¬ 


sion for 9 above to parameterise P in terms of 9, thus 
yielding 


/2 = 8 


1 — 2 tan ( - 


(56) 


P = 6V2 


Ni { 2 tt - n 0 ) 
^3 A/ 


(1 - A ) 


(57) 


After substituting these results into the Regge equa- 
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0 = 


Ni 


{27T -nO) A,; , 1/2_£/4_. 

2 64 


n li'^l 


N.A 


,2 


12^2 


8^2 3 - ii2 


3 ■■ 


A^2 II 
2V2 1 - 


11+6 [- 1-1 


(58) 


Once again, the first equation is the Hamiltonian con¬ 
straint while the second is the evolution equation. From 
the Hamiltonian constraint, we see that one need only 
specify l(t = 0) = Iq as initial data on some initial Cauchy 
surface, and then l{t = 0 ) will follow from the constraint. 

It can be shown that the Hamiltonian constraint is ac¬ 
tually a first integral of the evolution equation; a proof 
has been provided in Appendix [BJ This implies that the 
Hamiltonian constraint is sufficient to determine the evo¬ 
lution of I (t ), so we shall henceforth work with (1571) only. 

Using (|77)) . we can express I (t) and I (t) parametrically 
in terms of 6 to obtain 


2 _ A^i (27r - n9) 
N 3 A tan (i6)) 


(59) 


For the strut-length to be time-like, that is, for < 

0, we require 9 > tt/S, and for I to be real, we require 
9 < 2arctan(l/-\/2). Hence, 9 must lie in the range 


— < 9 <2 arctan 

O 



(60) 


Now that we have l{t) and /(t), we can determine R{t). 
Recall that we chose the integration constants in a{t) so 
that a = 0 when t = 0 . At t = 0, we have that = 3/A. 
On the other hand, dR/dt will be zero when I = 0, which 
happens when 9o = 2, arctan(l/-\/2). Inserting 9o into 
dini), we have that 


_ 3 2 V 2 N 1 
~ A Ns 


2 'k — 2 n arctan 



Therefore, we find that 


m 


/ Ns 

1 

\ 2 V 2 N 1 

27r — 2 n arctan 


1/2 

lit). 

( 61 ) 


C. Local variation of the parent models 

For comparison, we shall now derive the Regge equa¬ 
tions by locally varying the action. As mentioned earlier, 
this requires fully triangulating our skeleton; each trape¬ 
zoidal hinge will now be divided by a diagonal into two 
triangular time-like hinges. We shall label the lower tri¬ 
angular hinge by A and the upper by B, as depicted in 
Fig. El 


k+l 



Figure 5. A diagonal di divides the time-like hinge into a lower 
and upper triangular hinge, labelled A and B, respectively. 
The two struts on the sides of the hinge are now considered 
independent quantities. 


After triangulation, the skeleton’s hinges now comprise 
the time-like triangular hinges, the space-like triangu¬ 
lar faces of the tetrahedra, and triangular hinges formed 
by two diagonal edges for its sides and one tetrahedral 
edge for its base. This last set of hinges, an example 
being ABC', are new in that they have no counterpart 
in the global skeleton, unlike the first two sets. When 
the lengths of the diagonals are set to be equal, these 
hinges reduce to a set of isosceles triangles. It can be 
shown though that the corresponding deficit angles for 
these hinges are zero, so by virtue of the Schlafli identity, 
they will not contribute to the Regge equations. Thus 
we can drop them from the action. 

Therefore the Regge action for the triangulated skele¬ 
ton can be written as 

8^ Slocal = + E - A E ’ 

ie{4-blocks} 

(62) 

where, again, the superscripts A and B denote the lower 
and upper triangular subdivisions of the trapezoidal 
hinge as depicted in Fig. [2j Here, we continue using 
as a short-hand to denote the sum of the entire 4-block 
volume even though the 4-block is actually triangulated. 

We shall focus on obtaining the Hamiltonian constraint 
first, that is, on varying with respect to the strut-lengths. 
Locally varying with respect to an arbitrary strut-length 
rrij, we obtain an equation of the form 




drrij 


' i9m, 


drrii 


( 63 ) 


where 
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dK 


(4) 


druj 


= NaA'B'C'D‘ 


dVol 


AA'B'C'D' 


drrij 


- + NabB'C'D> 


dVolABB'C'D' 


drrij 


+ NaBCC'D‘ 


dVolABC C'D' 
dtrin 


-+Nabcdd' 


dVol^BCDD' 

dm 4 


and where Vol^ denotes the volume of 4-simplex X situ¬ 
ated between and Si+i. Na and Nb are the numbers 
of lower and upper time-like triangular hinges meeting 
at mj, while Naa'B'C'D', Nabb'CD' , Nabccd', and 
Nabcdd' are the numbers of 4-simplices corresponding, 
respectively, to AA'B'C'D', ABB'C'D', ABCC'D', and 
ABCDD' meeting at mj as well. This equation only 
depends on geometric quantities involving the time-like 
hinges and the volumes of the triangulated 4-block. 

We first consider the pair of triangular time-like hinges. 
Using Heron of Alexandria’s formula of classical antiq¬ 
uity, we can express their areas Af and Af in terms of 
the edge-lengths giving 



( 65 ) 


Differentiating with respect to their respective strut- 
lengths, we obtain 


dAt 

dmf 

dA^ 

drn^ 


A 

W 

m^ 

Slf 


+ df - 



( 66 ) 

(67) 


If we now substitute in the strut and diagonal lengths 
given by m and nil), and then take the continuum 
time limit, these derivatives simplify to 


dA^ _ dA^ 

dm^ dmf 



-l- 0{dt ), 


( 68 ) 


where again, in the context of continuum time, the over¬ 
dot denotes a time-derivative. 

After edge-lengths of the same type have been set 
equal, the deficit angles df' and i5f of the time-like trian¬ 
gular hinges become identical to the deficit angle of 
the original trapezoidal hinge; that is, = df = 
Setting the edge-lengths equal makes the two triangu¬ 
lar hinges be co-planar both with each other and with 


the original trapezoidal hinge; thus the 4-blocks meet¬ 
ing at the triangular hinges would be flat; the unit nor¬ 
mals of the triangulated faces would be identical to the 
unit normals of the original faces; and the dihedral an¬ 
gles between the triangulated faces would be identical to 
the dihedral angles between the original faces. Since the 
number of faces meeting at the triangulated hinge is the 
same as the number of faces at the original hinge, the 
dehcit angles for the triangulated hinge and the original 
hinge are identical. Additionally, since the dihedral an¬ 
gles are unchanged, then in the continuum time limit, I 
would still be given by (1551) . 

We next consider the volumes of the 4-simplices. To 
compute these volumes, we shall use the comparatively 
modern Cayley-Menger determinant instead, which gen¬ 
eralises Heron’s formula from areas of triangles to vol¬ 
umes of n-simplices. Suppose we label the vertices of 
an n-simplex with numbers 0 to n; then the simplex’s 
volume can be computed by the formula 




■(_l)n-si 

2"(n!)2 


1/2 


det(H) , 


(69) 


where H is a symmetric matrix given by 

fori,/>l, 

0 for z = j = 1, (70) 

1 otherwise, 


and where hj = Iji is the distance between vertices i and 
j. So for the volume of a 4-simplex, we have 

VolA) = ^^^-det{B), ( 71 ) 

and 



1 

1 

1 

1 

1 \ 

1 

0 

B 

‘01 

B 

‘02 

B 

‘03 

B 

‘04 

1 

‘oi 

0 

B 

‘l2 

B 

‘l3 

B 

‘14 

1 

‘‘02 

B 

‘l2 

0 

B 

‘23 

B 

‘24 

1 

B 

'os 

B 

‘l3 

B 

‘23 

0 

ll. 

Vi 

B 

‘04 

B 

‘l4 

B 

‘24 

B 

‘34 

0 / 


We now differentiate each 4-simplex volume with re¬ 
spect to its strut-length. In our 4-block, each strut will 
be an edge of exactly one of the four 4-simplices, and each 
of the four 4-simplices will be attached to exactly one of 
the four struts. When we differentiate each volume with 
respect to its associated strut-length, the resulting ex¬ 
pression simplifies greatly if we then take the continuum 
time limit; therefore we shall present only these contin¬ 
uum time expressions here. In the continuum time limit, 
we find that 
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dVol^A'B’C'D' 

dmf^' 


dVol^BB’C'D' 


dmf^' 

24^2 V 8 


dVol^BCC'D' 

' d^' 

1/2 

+ 0 {dt). 


dmf^' 


(73) 


With the continuum time form of the relevant geomet¬ 
ric quantities, we can now take the continuum time limit 
of the constraint equation (1631) to obtain a differential 
equation for I (t). Before doing this though, we can make 
some further simplifications to the continuum time form 
of (|6^ . We saw in the continuum time limit that all 
derivatives of the time-like hinges’ areas became iden¬ 
tical, as did all derivatives of the 4-simplices’ volumes. 
Using this knowledge, we can significantly simplify the 
continuum time form of (16311 to 


A AT ^ 

edges/vertex ^ .^^tetrahedra/vertex" 


The constants .hedges/vertex '^.nd A^tetrahedra/vertex ^^-re the 

numbers of edges and tetrahedra, respectively, meeting 
at any single vertex; they follow from the fact that each 
strut corresponds to the world-line of a vertex, each time¬ 
like hinge at that strut to an edge at the vertex, and each 
4-simplex at that strut to a tetrahedron at the vertex; 
therefore, the total number of triangular hinges meeting 
at the strut is identical to the number of edges at a vertex, 
-hedges/vertex j ^nd the total number of 4-simplices at the 
strut is identical to the number of tetrahedra at a vertex, 

-^tetrahedra/ vertex ■ 

Substituting all geometric quantities into the above 
equation, we obtain 




I 2 V 2 

A 


N, 


edges/vertex 


N, 


tetrahedra/vertex 


(27r — nO) 


However, Verges/vertex 'ind Ntetrahedrsi/verte'x 'ire related 
to Ni and N 3 by the relations 


-^edges / vertex 


2fVi 


and 


-^tetrahedra/vertex ; 

iVo 

where Nq is the number of vertices in a Cauchy surface 
and is also given by Table HI We therefore recover equa¬ 
tion (113; thus in this case, local variation yields the same 
Hamiltonian constraint as global variation. 

By locally varying action (l62ll with respect to k, we 
can also obtain an evolution equation identical to (1551) . 
To do this however, it turns out we must also make use of 
the Regge momentum constraints, that is, the equations 
obtained by locally varying the action with respect to 


each di- It can be shown that such equations are of the 
form 


o = aE 

3 


ddi 


(74) 


where the summation is over all 4-blocks containing the 
diagonal being varied. The area terms vanish from this 
equation, as it can be shown that dAf jddi = —dA^jddi 
for any triangulated trapezoidal hinge. 

To understand the relationship between the diagonal 
Regge equations and the evolution equations, let us first 
consider the variation of the 4-simplex volumes with re¬ 
spect to the diagonals and the tetrahedral edges. These 
derivatives actually depend on which tetrahedral edge 
or diagonal is being varied. Consider the triangulated 
4-blocks lying between Ei_i and Ei+i; if we vary the 
volumes with respect to the diagonals and then take the 
continuum time limit, we obtain 


ddf^' 

ddf^' 

ddf^' 

ddf^' 

ddf°' 

adf^' 


idt ,0 

24^2 

iSt ,0 

24^2 
iSt 
24^2 
iSt ,0 
24^2 
iSt ,0 
24^2 

iSt ,0 

24^2 


[\p)+o{df), 

[lp)+o{df), 

(i/2 + l)+0(rfi2) 


Clearly the diagonals give different derivatives. If we do 
a similar variation with respect to the tetrahedral edge- 
lengths, it can be shown that all six derivatives can be 
expressed in the form 


a(uE + u/1) 


dl 


xy 


iSt 



dd^y' ’ 


where is the length of the edge between vertices x and 
y and d^y is the length of the diagonal triangulating the 
world-sheet between Si and Si+i of edge . Thus each 
tetrahedral edge also gives a different derivative. 

Now if we vary the action with respect to an arbitrary 
edge , we obtain an evolution equation that can be 
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expressed in the form 

/ ,_ 5^(4) 

0=(RHSof(l58|)+A^^, (75) 


where the summation is over all 4-blocks containing di¬ 
agonal ; and by the diagonal Regge equation (TMl) . 
this summation is equal to zero. We thus recover the 
same evolution equation (155)) as we obtained from glob¬ 
ally varying the action. 

Since the diagonal equation is just a sum of different 
volume derivatives, we deduce from the form of the vol¬ 
ume derivatives above that the diagonal equation will 
necessarily have the form 


0 = -A 


iSt ,9 
24^2 



for some constants P and Q. Unfortunately, this equa¬ 
tion does not give any physically meaningful solutions 
as it implies that either I = 0 or P = — 8 Q/P. If P 
and Q have the same signs, then I would be imaginary. 
Even if P and Q have opposite signs, P would still be 
a constant, and relation (1551) for I would only equal this 
constant at a single value of d; this implies that there 
is only one moment in the evolution of the universe, as 
given parametrically by relations (I59p and (1561) . where 
such a diagonal equation could be satisfied. 

Furthermore, without knowing how the triangulated 
4-blocks fit together in the skeleton globally, we cannot 
determine the exact values for P and Q. The reason is as 
follows: a single diagonal will be shared by n 4-blocks; the 
diagonal might behave like a AB'-type diagonal in one 
4-block but like a AC-type diagonal in a neighbouring 
4-block, and these two 4-blocks will clearly have different 
contributions to the sum in dZH); hence it is even con¬ 
ceivable that different diagonals may give different Regge 
equations. 

We can however obtain one constraint equation by 
summing the diagonal Regge equations over all diagonals 
in the Cauchy surface Ej, where Ej denotes the Cauchy 
surface at time t in the continuum time model; it then 
follows that 


o=E 

in St 


E 

yW 

at 


dCv' 


= E 

y(4) 

in St 


E 

d^v' 
in V^^'> 


dd^y' 


= -n,AL,2 

24V2 



D. Relationship between the global and local 
Regge equations 

As Brewin has pointed out, the global Regge equa¬ 
tions can usually be related to the local Regge equations 
through a chain rule. A necessary condition is that the 
global and local Regge actions be identical. The summa¬ 
tions over hinges in the global action (l35ll and in the local 
action (15^ are identical for the following reasons: the di¬ 
agonal hinges do not contribute to the local action; the 
area of each trapezoidal hinge would always equal the 
areas of its two constituent triangular hinges while the 
deficit angles would always be identical; and both the ar¬ 
eas and the corresponding deficit angles of the triangular 
space-like hinges would always be identical in both the 
original and triangulated skeletons. The volume com¬ 
ponents of the two actions are also identical since the 
volume of a CW 4-block should equal the sum of the vol¬ 
umes of its four constituent 4-simplices. Therefore in this 
case, the global and local actions are indeed identical. 

Then by use of the chain rule, we can relate the global 
Regge equations to the local Regge equations. We begin 
with variation with respect to the struts. We can express 
the global variation of the common action S with respect 
to a global strut-length m as 

0 = — = V+ V —— ITBI 

dm ^ dm^ dm ^ ddi dm ’ 

i ^ i 

where the summations are constrained to the region be¬ 
tween a single pair of Cauchy surfaces Ej and Ej+i, and 
where m\ denotes the length of a local strut and di the 
length of a local diagonal. Since we shall be setting all 
strut-lengths equal, then dmj/dm = 1 for all i. Ad¬ 
ditionally, as we saw above, dS/dmf is 0(1) to leading 
order in dt. Thus the first summation has an overall lead¬ 
ing order of 0(1). On the other hand, it can be shown 
that 


ddi "tn 

dm d ’ 

which becomes m/l in the continuum time limit. How¬ 
ever m has a leading order of 0 {dt) in this limit, and 
as we saw above, dS/ddi has a leading order of 0{dt). 
Thus the second summation has an overall leading order 
of 0 [dP), which is higher than that of the first sum¬ 
mation. It therefore does not contribute to the Regge 
equation at leading order, and we can consequently sim¬ 
plify (176)) to just 


0 = 


dS 

dm 


E 


dS 

dmV 


(77) 


If Z ^ 0, then this implies that P = — |, which, as we 
have just remarked, is clearly non-physical. 

We shall discuss possible reasons for the local model’s 
unviability later on. 


clearly indicating that any solution of the local Regge 
equation will automatically be a solution of the global 
Regge equation as well. And as we saw above, the two 
solutions are in fact identical. 
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Using the chain rule, we can also express the global 
variation of S with respect to the tetrahedral edge- 
lengths I as 


^ ^ ^ 55 azf 

~ f)i ~ 




dS ddi 




dS dd,i 


5/f dl ^ dd, dl ^ dl 


(78) 

where the summations are constrained to a single Cauchy 
surface Y,j as well as the regions between and its neigh¬ 
bours Ej_i and Sj+i, and where denotes the length of 
a local tetrahedral edge. Because the tetrahedral edge- 
lengths will all be set equal, then dlj/dl = 1 for all i. Ad¬ 
ditionally, it can be shown that ddi/dl = ddi-i/dl = 1/2 
in the continuum time limit. Thus, we can simplify the 
above equation to 







(79) 


where we have made use of the fact that J^idS/ddi = 
'^^dS/ddi-i in the continuum time limit. However, it 
can be shown that in this case, both dS / dlf and dS/ddi 
are 0{dt) at leading order. Thus in contrast to the situa¬ 
tion with the struts, a solution to 0 = dS/dlf by itself is 
not sufficient to be a solution to 0 = dS/dl ; we must also 
satisfy 0 = dS/ddi. This is what we saw above, where 
we were able to recover the global evolution equation (1551) 
from its local counterpart (ESI) only when we also made 
use of the diagonal Regge equation (Tfil) . 


E. Initial value equation for the parent models 

In (3-1-1 )-formulations of general relativity, one would 
customarily determine the entire space-time by specify¬ 
ing a set of data on some initial Cauchy surface Eg and 
then evolving that data forwards in time to determine the 
rest of the space-time. Naturally, the evolution equation 
would be derived from the Einstein field equations. It has 
been shown that the required initial data consists of 
the first and second fundamental forms, h and x; the 
former corresponds to the projection of the metric g into 
Eg and effectively determines the 3-dimensional intrin¬ 
sic curvature of Eg; the latter effectively determines the 
extrinsic curvature of Eg within the overall space-time. 

There is a set of constraint equations that the initial 
data must satisfy to be consistent with the Einstein field 
equations. Let us express the Einstein field equation in 
the form 


G = SttT, 

where G is the Einstein tensor and T the stress-energy 
tensor. Let n denote a field of normalised one-forms ev¬ 
erywhere orthogonal to Eg. By making use of the Gauss 
equation 

^ ^ = RafijS ^ ^ ^ ^ p Xpo- Xi^p T Xpp X^'<y^ 

(80) 


which relates the 3-dimensional intrinsic curvature 
^^^Rfivap of Eg to its extrffisic curvature x E^nd its 4- 
dimensional intrinsic curvature Rap-yS, we can express 
the relation G (n, n) = SttT (n, n) as 

(3)77 + ix^'' X"" hpp K, = 167rp, (81) 

where ^^'>R is the 3-dimensional Ricci scalar of Eg and p 
is the energy density of the matter source as measured 
by an observer co-moving with respect to Eg. Equation 
(|8T|) gives the first constraint equation; it is actually the 
Hamiltonian constraint of the ADM formalism, where it 
is customarily derived by extremi sing the ADM action 
with respect to the lapse function [^. Let {ui}, for i = 
1,2,3, denote a set of normalised basis vectors tangent to 
Eg, and let | denote covariant differentiation with respect 
to the metric connection implied by h. By making use of 
the Gauss-Codazzi equation 

Rap n'^hP^ = x%|<, - X‘^a|p> ( 82 ) 

which relates the extrinsic curvature x of to its 4- 
dimensional intrinsic curvature in the form of the Ricci 
tensor we can express the relation G(n, u^) = 

87rT (n, u^) as 

Ku - x'^'^i^ u'^i = nP (83) 

which is actnally a set of three equations, one for each 
i. This gives the rest of the constraint equations; these 
are the momentum constraints of the ADM formalism, 
where they are customarily derived by extremising the 
ADM action with respect to the shift functions . 

Quite often, the initial surface Eg is chosen to be the 
snrface at a moment of time symmetry, that is, the mo¬ 
ment when the surface’s extrinsic curvature, as given 
by the second fundamental form Xi vanishes. In this 
case, the momentum constraints would vanish while the 
Hamiltonian constraint would simplify to 

(3)i? = IGttp; (84) 

this is known as the initial value equation at the moment 
of time symmetry. 

When there is a cosmological constant, the Einstein 
field equations take the form 

G -f Ag = 87rT. 

However we can always absorb the cosmological constant 
term Ag into T; it effectively acts as a perfect fluid source 
where pA = —PK = A/87r; then the initial value equation 
for a vacuum A-universe can be expressed as 

(3)i7 = 2A. (85) 

We shall now demonstrate that when this equation is 
applied to the time-symmetric Cauchy surface of our A- 
FLRW Regge model, the equation is satisfied. However, 
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the initial value equation (1841) and the Einstein field equa¬ 
tions, from which it is derived, will only be satisfied in 
an average manner on a Regge Cauchy surface. Curva¬ 
ture in the surface is concentrated at the hinges only, yet 
matter can be distributed away from the hinges where 
the skeleton is flat; thus the two sides of the equation 
will not agree in a point-wise manner. This contradiction 
arises because the Einstein field equations actually apply 
to smooth manifolds rather than Regge skeletons; they 
come about by varying the Einstein-Hilbert action when 
the underlying manifold is smooth rather than discrete. 
Thus by using the Einstein equations in this manner, we 
are effectively varying the Einstein-Hilbert action on a 
smooth manifold first and then applying the resulting 
field equations on a discrete manifold afterwards. The 
standard approach in Regge calculus is to use a discrete 
manifold from the very beginning, with the field equa¬ 
tions obtained being different as a result. Clearly, the 
two approaches are not equivalent. 

The parent models have a moment of time symmetry 
at the point of minimum expansion. This happens when 
1 = 0, corresponding to a dihedral angle of 

00 = 2 arctan f —| = arccos 

V V2j 

and from the Regge equation dSi), the edge-lengths 
would then be 



=6V2 ^j{2tt -nOo). 


( 86 ) 


which is exactly identical to (l86l) provided S = (27r —n0o)- 
Since n triangular faces do meet at any single edge, 
then the deficit angle 8 will have the form (27r — n0o), 
where 0o is the 3-dimensional dihedral angle between 
triangles in the Cauchy surface. Thus to complete our 
proof, we must show that 9q = 6q. Consider a typical 
tetrahedron in with vertices A, B,C, D, and assign 
the vertices to have co-ordinates identical to the spatial 
co-ordinates of their counterparts in (IT3l) . Edge AB has 
faces ABC and ABD meeting at it; the unit normal to 
ABC is 


f>‘ABC — (Oj 0, 1), 


and the unit normal to ABD is 


nABD 



2^2 l\ 


thus the dihedral angle between the two faces is given by 


cos 00 = i^ABC • il-ABD 
_ 1 

“ 3‘ 


We therefore see that 0o = 0o = arccos(l/3), and hence 
our models do satisfy the initial value equation at the 
moment of time symmetry. 


F. Discussion of the parent models 


We shall show that this specifically is consistent with the 
initial value equation. 

If a Regge Cauchy surface is sufficiently uniform such 
that there is a well-defined ‘volume per vertex’, then as 
Wheeler noted [l^, an average for can be given by 


(3)^ _ _ Si _ 

‘volume per vertex’ ’ 


(87) 


where the summation is over all edges radiating from a 
single vertex. In this 3-dimensional Regge Cauchy sur¬ 
face, the tetrahedral edges are now the hinges, and Si 
is the 3-dimensional deficit angle corresponding to edge 
li. As the parent CW Cauchy surfaces are clearly very 
uniform, there is a well-defined ‘volume per vertex’ given 

by 


‘volume per vertex’ 


Ns 

— (volume of one tetrahedron) 
No 


Ns 


6^2 iVo 


'■o- 


Since there are 2A^i/A^o edges radiating out from any 
single vertex in the parent Cauchy surfaces, the summa¬ 
tion in dUI) can be expressed as = ‘2{Ni/No)lo 8, 

where 8 is the common deficit angle of all edges. There¬ 
fore, the initial value equation ((85ll . when applied to our 
model, becomes 


Before examining the behaviour of our global Regge 
models, we shall first postulate on the reasons for the 
local models’ failure. After we have obtained the Regge 
equations from local variation, the standard approach in 
Regge calculus would be to specify a set of initial data 
on a single Cauchy surface and then determine the edge- 
lengths on all subsequent surfaces using the Regge equa¬ 
tions alone. Rather than doing just this, we have ad¬ 
ditionally constrained edges on each subsequent surface 
to be identical, in accordance with the CW geometric 
constraints. As a result, the only variable left for the 
Regge equations to determine is the overall scaling of the 
edge-length on each surface. The extra constraints were 
motivated by the Copernican principle, as we expected 
each Cauchy surface to be homogeneous and isotropic 
like the FLRW Cauchy surfaces they are intended to ap¬ 
proximate. The constraints implicitly assume that if we 
evolve from a Cauchy surface with identical edges, our 
subsequent surfaces will continue having identical edges 
because of the Copernican principle; we had assumed 
this would be the outcome even if we did not explicitly 
impose the constraints, so the constraints were really ex¬ 
pected only to simplify calculations. The Regge equa¬ 
tions presented in (1571) and (1751) were obtained after the 
constraints were imposed. 

In the global model, all Cauchy surfaces did possess 
Copernican symmetries as we could readily swap any pair 


I 


2 _ 
0 — 


6^2 


Ni 

Ns A 


S, 


( 88 ) 
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U{t) 


Figure 6. The rate of expansion of the universe’s volume 
dU/dt versus the volume U itself for the FLRW universe and 
the three different Regge models. The Regge universe volume 
is given by the sum of the volumes of the Cauchy surface’s 
constituent tetrahedra. 



Uit) 


Figure 7. The rate of expansion of 3-sphere volumes dU/dt 
versus the volume U itself for the FLRW universe and the 
three different Regge models. We have used R{t) as defined 
in ([5T|l to be the Regge universes’ 3-sphere radii. As R{t) 
was defined for all models to equal the FLRW scale factor 
a{t) when dR/dt = a = 0, then the volumes for all Regge 
models should equal the volume for the FLRW universe when 
dU/dt = 0; hence all graphs above coincide at dU/dt = 0. 

of vertices, tetrahedral edges, or struts on a surface with¬ 
out really changing the surface itself. However this was 
no longer the case in the local model when we triangu¬ 
lated the Cauchy surfaces, as the diagonals seemed to 
disrupt homogeneity. Not all vertices, for example, were 
connected to diagonals in the same way; thus they could 
not be swapped without non-trivially altering the surface 


itself: in a 4-block, like the one depicted in Fig. [T] and 
Fig. [21 vertices D and A' would not be assigned any di¬ 
agonals while A and D' would each be assigned three. 
Because we no longer have perfect homogeneity, our ex¬ 
pectation for edges to remain identical under evolution 
was no longer well founded. Rather, we should allow the 
initial surface to evolve according to the unconstrained 
Regge equations; without the additional constraints, we 
should expect to have a different equation for each edge 
on the surface. 

There may perhaps be a third method of varying the 
Regge skeleton, lying somewhere between a completely 
global variation and a completely local one, such that 
each edge of the original CW skeleton could be varied 
individually without having to break the symmetries in¬ 
herent in the Cauchy surfaces. We shall refer to this third 
approach as semi-local variation. When a single edge of a 
4-block is being varied, it is possible to impose constraints 
on the 4-block’s internal geometry, without needing to in¬ 
troduce extra independent edges, such that the geometry 
would still be well-defined under variation. For instance, 
one could constrain each block’s internal diagonals to be 
specific functions of the external edges; then when an ex¬ 
ternal edge gets varied explicitly, the diagonals would get 
varied implicitly in accordance with the constraints. The 
idea behind semi-local variation is to impose such con¬ 
straints on the 4-blocks attached to the edge being varied; 
after variation, one would then impose whatever further 
constraints on the skeleton that are necessary to make 
all other 4-blocks consistent with CW 4-blocks. We shall 
refer to the first set of constraints, the ones imposed on 
the 4-blocks around the varied edge, as local constraints 
and the remaining constraints as global constraints. The 
local constraints, however, must be chosen in a way that 
keeps all tetrahedral edges identical, otherwise the model 
may still not be viable. We leave to future consideration 
whether such a choice of constraints is possible. 

We now turn to comparing the global Regge models 
against the continuum model. We shall consider both 
the 3-sphere radii of the different models as well as the 
total volume of the universe. In Fig. [51 we have plotted 
the rate of expansion of the universe’s volume against 
the volume itself for each of the models. The volumes of 
Regge universes were given by dsni), while the volume of 
the FLRW universe was given by ©■ We see that the 
Regge models have produced the correct qualitative dy¬ 
namics of the space-time; for instance, they do not show a 
universe that eventually collapses back in on itself, nor do 
they show any instabilities in the evolution; rather, they 
show a universe that continues expanding indefinitely, 
like the continuum space-time being approximated. In 
terms of accuracy, we see that as the number of tetra¬ 
hedra increases, the model’s accuracy improves, with the 
600 tetrahedra model matching the FLRW space-time 
especially well when the universe is small. 

For comparison, we have also plotted analogous graphs 
in Fig. |7| where the Regge universes’ volumes were taken 
to be volumes of 3-spheres of radius R{t), as given by 
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Figure 8. The expansion rate of the universe’s radius versus the radius itself for the FLRW model and the three different Regge 
models, where the Regge radius is taken to be (a) R{t), as given by (I61II . (b) the average radius R{t), as given by (I28II . (c) the 
radius R{t) to vertices, as given by (1211) . (d) the radius Ri{t) to the centres of edges, as given by (1221) . (e) the radius i? 2 (t) to 
the centres of triangles, as given by (1241) . and (f) the radius R 3 (t) to tetrahedral centres, as given by (1261) . 
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St 


Figure 9. A schematic diagram, projected onto a 1+1 plane, 
of a 4-block with space-like struts. The dashed lines corre¬ 
spond to null curves that originate from the 4-block vertices 
in Ef. These lines separate out the region of the 4-block inte¬ 
rior that can be reached by causal curves from Ef, the white 
region, from those that cannot, the shaded region. Because of 
this causal structure, no past-directed causal curve from any 
vertex in T,t+dt will intercept Ef, and therefore, Ef+dt lies 
outside the future domain of dependence of Ef. For Brewin’s 
case, the above figure should be vertically inverted, since for 
a contracting universe, Ef+df would be smaller than Ef. 



Figure 10. The ratio of the Hubble radius I/Hq to the tetra¬ 
hedral edge-length I {t) versus the edge-length for each of the 
three Regge models; the graphs terminate when the struts 
turn null. The Regge Hubble parameter Hq is given by 


(ICTI) . In these graphs, the Regge models also very closely 
approximate the FLRW universe at low volumes before 
diverging as the universe gets larger; the approximation 
again improves as the number of tetrahedra increases. 

In Fig. [51 we have plotted the expansion rate of the 
universe’s radius against the radius itself, with each sub¬ 
figure using a different measure of 3-sphere radius for the 
Regge models. In all cases, the approximation to FLRW 
again improves as the number of tetrahedra increases but 
gradually diverges as the universe expands. The figures 
also reveal that radius R{t) gives the best approxima¬ 
tion to the FLRW model. This was somewhat expected 
given that R{t) was deliberately defined so that it would 
match a(t) exactly at the point of minimum expansion, 
the point which corresponds to the beginning of all the 
graphs. However, the radius R 2 {t) to the centres of tri¬ 
angles also gives a very good measure, as Fig. I5el reveals: 
this clearly indicates that centres of the triangles lie very 
close to 3-spheres of radius R, but it is not clear if there 
exists any underlying reason for this. 

All of our graphs terminated at an end-point on the 
right; this corresponds to the moment when the time-like 
struts turn null. Brewin remarked on a similar feature 
in his dust-filled FLRW models [l^: he obtained analo¬ 
gous graphs that also terminated at an end-point when 
the struts turned null. However, Brewin was consider¬ 
ing a closed universe, and the end-points appeared while 
the universe was contracting, whereas we are consider¬ 
ing an open universe, and the end-points appear while 
it is expanding. Brewin noted, in his case, that at the 
point when the struts turn null, the surface Ej would no 
longer lie in the past domain of dependence of the sur¬ 


face Tit+dt', there would be points in 'St+dt which cannot 
be reached by any causal curves from any point in E*. 
Since our Cauchy surfaces are expanding instead, it is 
the reverse that happens; that is, T,t+dt no longer lies in 
the future domain of dependence of Ej, as illustrated by 
Fig. |9l Brewin suspected, in his case, that the end-point 
signalled the local curvature had become too large for the 
approximation to handle. Nevertheless, in both Brewin’s 
models and ours, increasing the number of tetrahedra 
does postpone the appearance of the end-point. In his 
case, he was able to reach smaller volumes with a larger 
number of tetrahedra, while we are able to reach larger 
volumes, as Fig. [7] shows. 

To investigate whether there is any relationship be¬ 
tween the struts’ turning null and the ratio of the Hub¬ 
ble radius I/iLo to the tetrahedral edge-length /(t), we 
have plotted this ratio against l{t) in Fig.lTOl The Regge 
Hubble parameter Hq is defined to be The fig¬ 

ure shows that the Hubble radius always falls below the 
edge-length before the struts turn null; thus there does 
not appear to be any relationship between the Hubble ra¬ 
dius and the edge-length at the moment the struts turn 
null. 

In all of our graphs, we noticed that the models diverge 
increasingly from FLRW as the universe gets larger. We 
believe this is due to the finite resolution of our models 


In the FLRW universe, Hubble parameter Hq is defined to be 
a/a. We can define the Regge parameter analogously, using 3- 
sphere embedding radii rather than the FLRW scale factor a(f); 
however, regardless of which embedding radius we use, the quan¬ 
tity will always reduce to 
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trying to approximate an ever-expanding universe. Re¬ 
gardless of the universe’s size, the number of tetrahedra 
in any given approximation is always kept fixed; therefore 
the resolution will degrade as the universe gets larger. 
This is consistent with our observation that the graph 
diverges more slowly from FLRW as the number of tetra¬ 
hedra is increased, as the resolution of the models with a 
higher number of tetrahedra are able to ‘keep up’ longer 
with FLRW. This also suggests that the CW formalism 
would be more accurate when approximating closed cos¬ 
mologies that do not expand indefinitely but rather at¬ 
tain a finite maximum size before collapsing in on them¬ 
selves; in this case, the approximation’s resolution would 
be better controlled. Nevertheless, as we saw above, even 
for an infinitely expanding universe, the approximation 
was still able to yield qualitatively correct dynamics; this 
suggests the formalism should still be helpful in studying 
the dynamics of other space-times in general, and the in¬ 
formation obtained would be especially invaluable if the 
exact space-time solution is unknown. 

IV. CHILDREN MODELS OF CLOSED 
VACUUM A-FLRW UNIVERSES 

From any parent model, one can always generate a 
secondary model by triangulating each parent tetrahe¬ 
dron into a set of smaller tetrahedra. Brewin has devised 
a method that can subdivide any tetrahedra such that 
not only can it be applied to parent tetrahedra to gen¬ 
erate secondary models, but it can also be applied to 
each child tetrahedron afterwards to generate even finer- 
grained models. Thus, Brewin’s scheme can in principle 
be applied indefinitely. However, we shall only consider 
the first generation of children models. 

Under Brewin’s scheme, each parent tetrahedron gets 
divided into 12 children tetrahedra. Seven new vertices 
are introduced, six at the parent edges’ mid-points, the 
seventh at the tetrahedral centre. If the parent ver¬ 
tices are labelled A, B, C, D, then {XY) will denote 
the mid-point vertex between any pair of parent vertices 
X,Y G {A, B, C, D}, while (ABCD) will denote the cen¬ 
tral vertex. A partially subdivided parent tetrahedron is 
shown in Fig. [TT] 

The scheme also introduces three distinct types of 
edges. Edges connecting parent vertices to mid-point 
vertices, such as A{AB), have length Ui. As there are 
six parent edges in a parent tetrahedron, so there are 
12 length-Ui edges. Edges connecting mid-points to mid¬ 
points, such as {AB){AC), have length Vi. Each face 
has three of these edges, and as there are four faces per 
parent tetrahedron, so there are 12 of these edges per par¬ 
ent tetrahedron as well. Finally, edges connecting mid¬ 
points to {ABCD), such as {AB){ABCD) have length 
Pi. Each mid-point contributes one such edge and each 
parent edge contributes one mid-point, so there are six 
length-pi edges per parent tetrahedron. No edges connect 
parent vertices to (ABCD). 


D 



Figure 11. A partially subdivided parent tetrahedron. Three 
of the mid-point vertices, (AB), (AD), and [BD), as well as 
new edges of length Vi connecting them have been shown. The 
parent edges AB, AD, and BD, have now been subdivided 
into two edges of length Ui each: A{AB) and {AB)B, A{AD) 
and {AD)D, and B{BD) and {BD)D. Also depicted is the 
central vertex {ABCD). 

As a result, there will also be three distinct types of 
subdivided tetrahedra, each with four members per par¬ 
ent tetrahedron, thus giving the total of 12 children tetra¬ 
hedra. Each Type I tetrahedron consists of an equilateral 
base formed by three mid-point vertices and an apex at 
a parent vertex; an example would be A{AB){AC){AD). 
There are four of these tetrahedra in a parent tetrahe¬ 
dron, one per parent vertex. Type II tetrahedra share 
the same equilateral base as Type I but have their apexes 
at {ABCD) instead. Because they share the same base, 
there is a 1-1 correspondence between Type I and Type II 
tetrahedra, so there are four Type II tetrahedra per par¬ 
ent as well. Finally, a Type III tetrahedron consists of an 
apex at {ABCD) as well and an equilateral base formed 
by the three mid-point vertices on a single parent face; an 
example would be {AB){AC){BC){ABCD). Since each 
parent face is associated with one Type III tetrahedron, 
there are four Type III tetrahedra per parent tetrahe¬ 
dron. We note that in terms of their edge-lengths. Type 
II and Type III tetrahedra are identical. The different 
vertices, edges, and tetrahedra in the subdivided parent 
tetrahedron have been summarised in Table imi 

Finally, each vertex’s world-line will generate a strut 
connecting one Cauchy surface to the next. As there 
are three sets of vertices, there will also be three sets of 
struts, with all struts in the same set sharing the same 
length. Following Brewin, we shall restrict our consider¬ 
ation to only models where all three sets are constrained 
to have the same length. 

Since all Cauchy surfaces in continuum FLRW space- 
time are identical to each other apart from a time- 
dependent scale factor a{t), by analogy, we shall require 
the subdivided Cauchy surfaces to be identical to each 
other as well apart from an overall time-dependent scale 
factor; that is, we shall require 

— = a and — = /3 Vi, 

Vi Vi 


( 89 ) 
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Table III. The vertices, edges, and tetrahedra of a subdivided 
parent tetrahedron. 




Number 

Simplex type 

Example 

per parent 
tetrahedron 

Parent vertices 

A, B, C, D 

4 

Mid-point vertices 

{AB), {AC), {AD), 
{BC), {BD), {CD) 

6 

Central vertex 

{ABCD) 

1 

Edge Ui 

A{AB) 

12 

Edge Vi 

{AB){AC) 

12 

Edge Pi 

{AB){ABCD) 

6 

Type I tetrahedra 

A{AB){AC){AD) 

4 

Type II tetrahedra 

{AB){AC){AD){ABCD) 

4 

Type III tetrahedra 

{AB){AC){BC){ABCD) 

4 


for some constants a and (3 independent of the Cauchy 
surface. This requirement is natural because our CW 
Cauchy surfaces are intended to approximate FLRW 
Cauchy surfaces; so if the only difference between two 
FLRW surfaces is an overall scaling gq Agq, then the 
only difference between the two CW surfaces approxi¬ 
mating them should be a re-scaling of all lengths by A as 
well. This is assuming that the two CW surfaces trian¬ 
gulate their respective FLRW surfaces in the same way; 
that is, no extra vertices, edges, or tetrahedra appear in 
one Regge surface but not the other. 

A. The 3-sphere embedding of children Cauchy 
surfaces 

The subdivided Cauchy surface will have a slightly dif¬ 
ferent embedding from that of the parent Cauchy surface. 
The main difference is that each of the three sets of ver¬ 
tices can lie on its own 3-sphere, as the three sets are inde¬ 
pendent of each other. Nonetheless, the three 3-spheres 
will share a common centre in E^. As the Cauchy surface 
expands or contracts, the radii of the three 3-spheres will 
increase or decrease correspondingly. Each 3-sphere can 
be parametrised using polar co-ordinates as given by HUD 
though with different radii. We shall denote the three 
radii by for the parent vertices, for the mid¬ 
points, and for the central vertices. As mentioned 
in the Introduction, our approach to this embedding will 
be completely different from that of Brewin [l^ , who in¬ 
stead constrained all three sets of vertices to lie on the 
same 3-sphere. 

We shall first embed a representative subset of the mid¬ 
point vertices. As with the parent model, we can always 
choose our co-ordinates such that one of these vertices lies 
at (0, 0,0). This vertex will have n + n nearest mid-point 


Parent vertex 



Figure 12. The upper and lower spokes of mid-point vertices, 
shown in solid lines, surrounding mid-point vertex 0 for the 
case of n = 4. One dimension has been suppressed. All 
vertices in the two spokes are situated on a sphere of radius 
Vi. Each m, m' pair of vertices are also nearest neighbours to 
each other, being located on the same parent triangle; they 
are hence separated by Vi. Each vertex is also separated from 
its two neighbours in the same spoke by a distance of Vi. 
Also depicted are the two nearest parent vertices and two 
of the nearest central vertices. The upper parent vertex is 
equidistant to each vertex of the upper spoke and to vertex 
0, while the lower parent vertex is equidistant to each vertex 
of the lower spoke and to vertex 0. The left central vertex 
is equidistant to vertices 0,3,4, 3', 4', and the right central 
vertex to 0,1, 2, 1', 2'. 

neighbours to it, where n is given by the final column of 
TablelH To help understand the positions of these nearest 
neighbours, we shall refer to Fig. 1111 Suppose our ver¬ 
tex at (0, 0,0) corresponds to the mid-point vertex {AB) 
in Fig. 111! then on each parent triangle containing this 
vertex, there will be two nearest mid-point neighbours, 
for instance {AD) and {BD) for triangle ABD in the fig¬ 
ure. Since n parent triangles share a parent edge, there 
will be n parent triangles containing vertex {AB)] thus 
we shall have a total of n nearest neighbours from {AD) 
and its analogues and another n nearest neighbours from 
{BD) and its analogues. This gives a combined set of 
n + n nearest mid-point neighbours. We shall refer to 
one of the two sets as the ‘upper’ set and the other as 
the ‘lower’ set. As these neighbours are all located at 
a common distance of Vi from the first vertex, they are 
situated on a common 2-sphere centred on this vertex. 
Hence, we can choose the x co-ordinate of these neigh¬ 
bours to be determined by just the radius of this 2-sphere. 
We denote this co-ordinate by x = Ao- 
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The upper and lower sets each form a spoke around 
the first vertex. An example of such spokes is illustrated 
in Fig. [12] for the case of n = 4. We shall label the upper 
vertices by 1 ... n and the lower ones by 1 '... n' such 
that if m G [1, n] is an (AZ?)-type vertex then m' G [1', n'\ 
would be the (i?Zl)-type counterpart located on the same 
parent triangle. With this choice of labelling, we note 
that each m and m! pair are nearest neighbours as well. 
We can choose the 0 co-ordinate to be such that the upper 
vertices are all located at 0 = 9q. Then if we choose 
vertex 1 to be at ^ = 0 , upper vertex m will be located 
at 0 = 2(m— l) 7 r/n. Having fixed the upper vertices, the 
lower vertices are constrained to be at d = tt — 0 o with 
vertex m! having the same (j) co-ordinate as its upper 
counterpart m. These results are summarised in Table 

CYl 


Table IV. The polar co-ordinates of a mid-point vertex and 
its 2 n nearest mid-point neighbours. 


Vertex 

X 

e 


0 

0 

0 

0 

1 

Xo 

00 

0 

2 

Xo 

00 

27r 

n 

n 

Xo 

00 

2(n — l)7r 
n 

T 

Xo 

1 

o 

0 

2 ' 

Xo 

. 1 

o 

27r 

n 

n’ 

Xo 

1 .. 

o 

2(n—l)7r 
n 


:= 


v^ 



l-cos^ ) 


(95) 


We see that once again, the ratio between the edge-length 
Vi and the 3-sphere radius R\ is independent of the 
Cauchy surface label i and hence of time. We see as 
well that the angular co-ordinates of the vertices are also 
independent of i, and therefore, as the edge-length Vi 
expands or contracts, the 3-sphere simply expands and 
contracts about its centre. 

We now consider the embedding of the nearest parent 
vertices. There will be two nearest parent neighbours to 
vertex 0. One of them will also be a nearest neighbour 
of all n upper vertices and the other of all n lower ver¬ 
tices, as illustrated in Fig. [T^ By symmetry, the upper 
parent vertex will be located on the 2 -dimensional plane 
containing vertex 0 but orthogonal to the plane contain¬ 
ing the upper n vertices. We can obtain one vector 
in this plane by taking the average of the n vertices and 
subtracting the position of vertex 0. This gives the vector 

(pW)'' = (-Z(2),l,o,o). (96) 


A second basis for the plane can be obtained by re¬ 
quiring orthogonality to both and any vector joining 
any two of the upper n vertices. This gives 


In terms of these two vectors, the location of the 
parent vertex is then 


f2) 

The parameters xoi So, can be related to the edge- 
length Vi- As mentioned above, the xo co-ordinate is 
fixed by the distance between vertex 0 and its nearest 
neighbours; this yields the relation 

u2=2(Af))"(l-cosxo). (90) 

From the distance between vertices 1 and 2, we have the 
relation 

v1 = 2 sin^ Xo sin^ 6 >o ^1 - cos , (91) 

and from the distance between 1 and 1 ', we have the 
relation 


= 4 



sin^ 6*0 cos^ xo- 


Solving these equations yields 


cos Xo = 


cos 6*0 = 


1 + cos f 
2 (l-cosf)’ 
/l_cos^y^^ 
U-cos^j 


(92) 


(93) 

(94) 


(1,0, 0,0) + A(i) ^ , 

since both vectors p*^^^ and are given with respect to 
the position of vertex 0. We also require be equidis¬ 
tant to vertex 0 and all upper n vertices; this constraint 
fixes to be 


A(i) 


1 + (Z(2))2’ 


(98) 


while remains a free parameter. 

It can be shown that q^^^ is actually parallel to the 
radial vector connecting the parent vertex to the centre 
of the 3-sphere, that is. 




(99) 


where q*^^^ is simply q^^^ normalised. We can therefore 
express more simply as 



R 


( 1 ) 


VT+W^‘- 


=(i,z(2),o,o), 


( 100 ) 


and radius 



is related to parameter by 



l-h(Z(2))2 



( 101 ) 
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Since the free parameter only changes along the 
direction of changing will only change the ra¬ 
dius of the 3-sphere. 

We can deduce the embedding of the nearest central 
vertices in a similar manner. As there are n parent tetra- 
hedra hinging on a parent edge, there will be a total of n 
central vertices that are nearest neighbours to vertex 0. 
Each of these n vertices will be located along the central 
axis to each vertex quadruplet of the form m, m', (m-|-l), 
{m+iy for all m G [1, n], with vertices (n-l-1) and (n-|-l)' 
being identified with vertices 1 and lb Two examples of 
central vertices are shown in Fig. [T^l We shall focus on 
the central vertex equidistant to vertices 0, 1, 1', 2, 2 '. 
Similar to the situation with the parent vertex, the cen- 

I 


tral vertex will be located on the 2-dimensional plane 
containing vertex 0 but orthogonal to the plane contain¬ 
ing vertices 1, 1', 2, 2'. A vector in this plane can 
be found by taking the average of vertices 1, 1', 2, 2' and 
then subtracting the position of vertex 0; this gives 




2-(Z(^))^ 


3-(Z(2))V )' 


. . 

The second basis can be obtained by requiring or¬ 
thogonality with and with any vector connecting any 
of vertices 1, 1', 2, 2'; this yields 




/ 2 - 

V3- (Z(2))^ 


1/2 




(103) 


Then in terms of these two vectors, the position of 
the vertex is 

= Af)(1,0, 0,0) + A(2) ^(2) (^(2)^'^. 

Requiring to be equidistant to vertices 0, 1, 1', 2, 2' 


yields 


(104) 

while remains a free parameter as well. 

As with the parent vertex, it can be shown that q^^^ 
is actually parallel to the radial vector pointing from the 
centre of the 3-sphere to the central vertex, that is. 


r(3) ^ i?(3)^(2)^ (105) 

where q*-^^ is simply q^^^ normalised. Therefore, we can 
also express as 



V 2 



( 2 -{Z^Y" \ 

V3-(^(^^)V ’^3-(Z(2))V ’ 


(106) 


where 


is related to 


by 


The relationship between Ui and is 


(107) 

Since the free parameter only changes along the 
direction of q'-^^ changing will only change the ra¬ 
dius of the 3-sphere. 

We have previously related Vi to in (1951) . We can 
also relate edge-lengths Ui and pi to the 3-sphere radii. 


w? = 


(«<■')’ 


?(3) 


(108) 


while that between pi and R\ ' is 

pI = {Rf)'-2RfRf^^2 - (Z(2))2+(i?f))^ (109) 

Since and are free parameters that determine 
the 3-sphere radii R^P and respectively, we can re- 
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express (|101l) and (11071) in the form 

R^^=aR^\ 


( 110 ) 


where the scaling factors d > 0 and j3 > Q now become 
the free parameters. This effectively means that and 
can be freely chosen for some initial Cauchy surface, 
and from (jl08() and (jlOOp . this choice would effectively 
determine Ui and pi. Therefore, the freedom to choose 
a and P is equivalent to a freedom to choose a and j3 in 
(l89l) for some initial Cauchy surface. 

By altering a or j3, we would expand or contract one 
of the 3-spheres relative to the others. The vertices on 
the altered 3-sphere would simply shift radially inwards 
or outwards, but not angularly. Because the shift is 
purely radial, a vertex on this 3-sphere would still re¬ 
main equidistant to its nearest neighbours on the same 
3-sphere, although that distance would change. Simi¬ 
larly, the vertex would remain equidistant to its nearest 
f2) 

neighbours on the R] 3-sphere but with the distance 
altered as well. 

As mentioned at the start of this section, all FLRW 3- 
spheres are identical to each other apart from an overall 
scaling a(t), and we approximate this symmetry by re¬ 
quiring ratios a and j3 as defined in (1891) to be constant; 
then the evolution of two sets of the Cauchy surface edge- 
lengths can be determined by the third set alone, which 
we shall take to be Vi. Since a and /3 are equivalent to a 
and /3, respectively, we can equivalently require that all 
of our Regge 3-spheres be identical to each other apart 
from an overall scaling; we would freely specify a and 
/3 for some initial Cauchy surface, but our requirement 
would constrain a and /3 to be the same for all subsequent 
surfaces. As a result, the evolution of two of our 3-sphere 

radii can be determined by the third radius alone, and 
f 2) 

we shall choose R\ to be that sole dynamical radius. 
The analogy between CW Cauchy surfaces and FLRW 
Cauchy surfaces is much clearer when working with these 
embedding 3-spheres and 3-sphere radii rather than with 
tetrahedral edge-lengths. However regardless of whether 
we work with radii or edge-lengths, there is only one dy¬ 
namical length parametrising the entire system. 


B. Child 4-block co-ordinates 

To facilitate the calculation of geometric quantities in 
the children models, we shall introduce a co-ordinate sys¬ 
tem similar to m- Our approach is a modification 
of Collins and Williams’ original approach and differs 
from that followed by Brewin, who uses non-Cartesian 
co-ordinates. Let us consider a typical tetrahedron of 
the model in Cauchy surface It will always have an 
equilateral base of edge-length Vi regardless of the tetra¬ 
hedron’s type. We label the base’s vertices by A, B, C 
and the apex by D. The three edges meeting at the apex 


will all be of identical length, either Ui or pi, and without 
loss of generality, we shall work with Ui. Then we have 
the freedom to set the co-ordinates for A, B, C, D to bcEl 


A = 
B = 






c= (o,^,o, 

D = (0,0, hi, di), 


( 111 ) 


where hi is the height of the tetrahedron and is given by 



This tetrahedron evolves to another in surface Si+i 
with vertices labelled A', B', C, D'. The vertices’ co¬ 
ordinates are determined by constraints on the edge- 
lengths and the requirement that the 4-block have no 
twist or shear. The edge-length constraints are that the 
lengths of struts AA', BB', CC be equal, the lengths of 
edges A'B', A’C, B'C be and the lengths of edges 
A'D', B'D', C'D' be Ui+i. The requirement of no twist 
or shear was also imposed on the parent 4-block; thus, 
the parent and child 4-block geometries are constrained 
in similar ways. Once again, the constraint on the strut- 
lengths and the requirement of no twist or shear can be 
considered analogous to a choice of shift and lapse func¬ 
tions in the ADM formalism. 

Evolution should preserve the equilateral symmetry of 
the tetrahedron’s base; therefore, the base should simply 
expand or contract uniformly about its centre. However, 
it may also undergo a vertical displacement 5zi, which is 
determined by the struts’ lengths. 

The requirement that apex D' be located at distance Ui 
from each base vertex is equivalent to the two constraints 
that D' lie at distance from the base’s centre and 
that the central axis connecting D' to the base’s centre 
lie orthogonally to the base. As base A'B'C defines a 
2-dimensional plane in a (3-|-l)-dimensional Minkowski 
space-time, the subspace orthogonal to it would be a 
(l-l-l)-dimensional plane, and the tetrahedron’s central 
axis can be oriented along any direction in this plane. 
Combined with the hrst constraint, the second constraint 
implies that D' will lie on a hyperbola in this (1-1-1)- 
dimensional plane; exactly where on the hyperbola it lies 
depends on the length of strut DD'. Therefore the axis of 


® We note tliat because of the underlying 3-sphere geometry, the 
apex of the neighbouring tetrahedron sharing base ABC will in 
general be Lorentz-boosted relative to vertices A, B, C, D in co¬ 
ordinate system lllllll ; its co-ordinates would in general take the 
form ^0, 0, —hi cosh Oi, iti -|- ihi sinh Qi^ , where hi is the height 
of the neighbouring tetrahedron and ai is a boost parameter. 
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the upper tetrahedron may in general be Lorentz-boosted 
relative to that of the lower tetrahedron. 

Thus the upper vertices’ co-ordinates are given most 
generally by 


_ I Vi+i V^+l , 


B' = 


2 ’ 2^3 

f Vi+l Vi+i 


V 2 ’ 2^3^ 


8zi , I, 


c = rfi+ij, 

D' = (0, 0, hij^i cosh'01 -I- 5zi, t0+i -I- sinh0i), 


where 0i is the relative boost between the upper and 
lower tetrahedral axes. 

However if the 4-block is to have no twist or shear, 
then we also require that the world-sheet generated by 
each tetrahedral edge between and be flat; in 

other words, the four vectors parallel to the four sides of 
this world-sheet must be co-planar. When this require¬ 
ment is imposed on any world-sheet involving D', such 
as ADA'D', and when the scaling relations (l8^ are im¬ 
posed, it can be shown that 0^ must be zero. Therefore, 
the upper tetrahedron’s co-ordinates can be simplified to 


A' = 
B' = 



Vz+l 

2^3 


, Sziy j , 


Vj+l 

2 


yj+i 

2V3 


, Szi , 


c = ^0, 

D' = (0,0, hi+i + 5zi, iti+i). 


(113) 


We can now deduce 5zi from the requirement that 
struts AA' and DD' have equal length; this gives the 
equation 


1 


6v^ + 5zf 



2 


5vi + Szi , 


where 6vi denotes 6vi := Vi+i — Vi. Solving this yields 


6zi = 




(114) 


C. Geometric quantities 

We shall now derive the geometric quantities relevant 
to the Regge equations. However, we shall only present 
those quantities relevant to the Hamiltonian constraint 
equation, that is, to varying the action with respect to 
the strut-lengths, as our plan is to use this constraint to 
study the behaviour of the children models, much like 
how we studied the behaviour of the parent models from 


the parent Hamiltonian constraint. Moreover, we shall 
only present the local variation of those quantities, since 
global variation can always be related to local variation 
by the chain rule. 

Before proceeding, we shall first remark on an impli¬ 
cation that the strut-length constraint has on the rela¬ 
tionship between a and (3. Using co-ordinates ()111|) and 
(11131) . we can express the length of a parent vertex’s strut 
as 


= 4 (a2 _ 1) 

where again, SU denotes SU := 0+i — ti. By simply 
swapping a for 0 in the above expression, we obtain an 
equivalent expression for the length of a central vertex’s 
strut, that is. 


mf = 


04 


4(0^-^) 


Svf — Stf. 


As we require all strut-lengths to be equal between any 
pair of consecutive vertices, it follows that 



I) 


04 

(0^ - ’ 


which implies that 


0 = a or 



(116) 


(117) 


This relationship between a and 0 implies that one of the 
three edge-lengths, and hence one of the three 3-sphere 
radii, can no longer be independent of the other two edge- 
lengths and radii. 

We shall now turn to presenting all geometric quan¬ 
tities relevant to the children models’ Hamiltonian con¬ 
straint. As with the parent models, the only relevant 
geometric quantities are the varied time-like hinges gen¬ 
erated by the tetrahedral edges’ world-sheets, the corre¬ 
sponding deficit angles, and the 4-blocks’ varied volumes. 

We begin with the variation of the time-like hinges. As 
with the parent models, the time-like hinges can again be 
triangulated in the same manner as depicted in Fig. [S] 
and we can again use equations (IMl) to (1571) to obtain the 
variation of the resulting triangular hinges with respect 
to their strut-lengths rrii. For the hinges generated by an 
edge of length u,, we find that 


^ a nii Vi+i + 

dmf- dmf 2SU _4’ 


(118) 


where Vi denotes Vi := SvijSti. For the hinges generated 
by an edge of length pi, we can swap a for /3 to obtain 


_ dAf^ _ I3mi Vj+i+Vj 

dmt dmf 25U /i^^2 _4' 


( 119 ) 
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Finally for the hinges generated by Vi, we find that given by 


_ rui {vj+i + Vj) 

dmf^ dmf 




( 120 ) 


As we shall ultimately take the continuum time limit 
of the Regge equations, we shall at this stage present the 
continuum time limit of the quantities above for later use. 
As 5ti —>■ 0, equation (11181) becomes 

av \ -—- I ; (121) 


dmf dmf 
equation (11191) becomes 


_ 4 

3 a2-i 


8/l“» a/lf’ „ / - 1 


1/2 


dinf i9mf 
and equation (I120D becomes 


^*2-4 

3 


( 122 ) 


dA^^ _ dAf^ 

dmf" dmf 


1 a „-,2 
4 a2- 


— 1 


1/2 




■ (123) 


COS 


^ = fii • n 2 


1 + 


4(a=-l) 


■.2 

4(a2-i) 4 


1/2 


3(4a2-l)- 


4(“"-3) 


1/2- 


(126) 

Let us next consider the dihedral angle at hinge 
ADA'D', which is generated by length-Mi edges. This 
hinge is also co-planar, and so the two triangular hinges 
that subdivide it will also have the same dihedral an¬ 
gle. The faces meeting at this hinge are ABDA'B'D' 
and AC DA'CD'. Face ABDA'B'D' repeats from be¬ 
fore and so has unit normal 112 , while face ACDA'CD' 
has unit normal 


/ 


3 (a2 - i) 


1/2 


-V3,l, 


- 1)’ 2P3(a2-i)y 

(127) 


We next turn to determining the deficit angles. In the 
4-block described by (11111) and (I113|) . there will be two 
different dihedral angles, one at hinges generated by the 
edges of length Vi and the other at hinges generated by 
the edges of length u^. 

Let us first consider the dihedral angle at trapezoidal 
hinge ABA'B', which is generated by a length-edge. 
Because this hinge is co-planar, the two triangular hinges 
that subdivide it will have the same dihedral angle as 
that of the original hinge. The two faces in the 4-block 
meeting at ABA'B' are ABCA'B'C and ABDA'B'D', 
and their unit normals are, respectively. 


= 1 - 


and 




4 (a2 _ 1 ) 


0 , 0 , 1 ,-L - , ^ Vi 


2x1 C - i 


(124) 


Q “ - 


3 (4a^ - 1) - 


4 -2 
a vf 


1 - 1/2 




X I 0,-2J3la^-i|,l,t 




(125) 


2x1 cP- 


Therefore the dihedral angle between the two faces is 


Thus the dihedral angle is given by 
AT) 

COS 0) ' = n2 • na 


3 ( 2 a 2 _ 1 ) + 
3(4a^-l)-77^ 


(128) 


There will be another pair of dihedral angles, anal¬ 
ogous to the two above, in the 4-block with edges of 
length Pi instead of U). These angles can immediately 
be obtained by swapping a for /3 in the two expressions 
above, thus giving 


cos6(p^ = 


and 


1-h 


4(/3=-p 




i_(^!zi)!.2 
4 4(/32-i'l^* 


1/2 r 


3(4/32-1)- 


-|l/2 


4(/3^-P 


h24=-i) + i^ 


cos6(p^ = 


3 ( 4/32 - 1) - 


(129) 


(130) 


4(/3^-p 


Therefore the child model has a total of four distinct 
dihedral angles. 

With these dihedral angles, we can now deduce the 
deficit angles at each hinge. We begin with hinges gen¬ 
erated by length-Mi edges. There will be n space-like 
triangular faces meeting at such an edge, with n again 
given by the last column of Table HI and hence there will 
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be n faces meeting at the hinge generated by that edge. 

Between each adjacent pair of faces is a dihedral angle of 
( 2 ) 

0l , and hence the hinge’s deficit angle is 

) = 27r - ^ (131) 

We next consider the deficit angle at hinges generated 
by length-rii edges. These edges are attached to all three 


types of tetrahedra, and each type contributes a different 
dihedral angle to the corresponding hinges’ overall deficit 
angle; each Type 1 tetrahedron contributes while 
each Type 11 and Type 111 tetrahedron contributes 0) . 
The number of Type I tetrahedra at a length-Wi edge is 
given by 


(No. of Type 1)(3 n^-edges per Type I) 

(No. of Ui-edges) 

(4 Type I per parent) (No. of parents) (3 z;i-edges per Type I) 
(3 Ui-edges per parent triangle) (No. of parent triangles) 

= 2 , 


r 


where ‘parent’ refers to a parent tetrahedron, and where 
we have used the fact that, regardless of the model, 


(No. of parent triangles) 
(No. of parents) 


= 2 ; 


(132) 


this last identity follows from the fact that there are al¬ 
ways four triangles per tetrahedron, but every triangle is 
always shared by two tetrahedra, so there are only half as 
many triangles per tetrahedron in a skeleton. By similar 
combinatorics, we can show that the number of Type 11 
tetrahedra is given by 




Type II 


(No. of Type 11) (3 Wi-edges per Type II) 


(No. of Ui-edges) 


= 2 , 


and the number of Type III by 

(No. of Type III) (3 n^-edges per Type III) 




Type III 


(No. of fi-edges) 


= 2 . 


Therefore the deficit angle is 

5 ( 2 ) = 271 - 20 ?^ - 


(133) 


Finally, we consider the deficit angle at hinges gener¬ 
ated by length-pi edges. Only Type II and III tetrahedra 
have such edges, and each will contribute a dihedral an¬ 
gle of . The number of Type II and III tetrahedra at 
an edge is 




Type II & III 


(No. of Type II & III per parent) (3 pi-edges per Type II or III) 


(No. of Pi-edges per parent) 


= 4, 


and therefore the deficit angle is 

=2tt- 40 ^^ 


time expressions again. For the 4-block generated by a 
Type I tetrahedron, the derivatives of all four 4-simplices’ 
(134) volumes simplify to 


The final geometric quantities to derive are the 4- 
simplices’ varied volumes. A 4-block generated by a 
tetrahedron can again be triangulated in the same man¬ 
ner as shown in Fig. [2l and the volumes of the four 4- 
simplices generated can again be calculated in the same 
manner as with the parent models, using equations (1691) 
to (17^ . Again, we vary each of these volumes with re¬ 
spect to their associated strut-length. The resulting ex¬ 
pressions also simplify greatly after the continuum time 
limit is taken, so we therefore present only the continuum 


dv^ 

drrii 




16V3 




+ 0{dt). 


(135) 

The 4-block generated by Type II and Type III tetrahe¬ 
dra are identical, and the derivatives of all 4-simplices’ 
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volumes simplify to 


fc III 

drrii 






1)2 ) + 0 (dt). 

(136) 


D. Varying the Regge action 

We shall now derive the Regge equations for the chil¬ 
dren models. We initially attempted to locally vary the 
Regge action with respect to the strut lengths, however 
this resulted in a mutually inconsistent set of equations. 
Because each child model has three distinct sets of struts, 
local variation would yield a total of three distinct con¬ 
straint equations, one for each set. Any of the constraint 
equations could then be used to determine the behaviour 
of our dynamical variable Vi, but with three constraint 
equations determining a single Vi, the system of equations 
risked being over-determined. This indeed happened: we 
obtained three constraint equations that gave mutually 
inconsistent relations for Vi regardless of which choice we 
made for /? in (11171) . Therefore, we can perform only a 
global variation of the action, and this is what we shall 
now present. 

The Regge action for a child model can be expressed 
as 


We can use the chain rule to express any globally var¬ 
ied quantity as a sum of locally varied quantities. For a 
common strut-length mi, the chain rule takes the form 


d dm^ d dmy d 9mf^ d 

_ ^ J __|_ _ J __j_ _ J _ 

drrii drrii drrii drrii 

d^_^ d^_^ ddf d 

drrii dim ^ dm^ ' 


Since all struts have equal length, then 


dmy 

drrii 


= 5, 


V/c. 


We also have that 


dd\ 


{k) 


drrii 


rrii 

j(fc) 


Sij = 0{dt) 


\/k. 


However we found that the leading order of the Regge 
equation 0 = dS/dmy was 0 ( 1 ) for all k, so once again, 
the diagonal derivatives do not contribute. Hence to lead¬ 
ing order, the chain rule can be simplified to the form 

d d d d 

drrii ~ dm^ dm^ dinf ' 



E[ 


iVArea 1 (^Al ^(1) 

+ iVArea 2 ^B 2 ^ ^( 2 ) 

+ iVArea 3 ^B3^ ^0) 


- 4N3A{vy+2V, 


.(4) 


"z.II & III 


(137) 


We have used this chain rule to globally vary the Regge 
action with respect to nii to obtain the Hamiltonian con¬ 
straint for the children models, and we have then taken 
the continuum time limit. The resulting equation is 


4-\/3 




-1 


where the summation i is over the Cauchy surfaces, V 3 is 
the number of parent tetrahedra, and is the volume 
of the entire 4-block of a Type X tetrahedron, that is the 
combined volume of all four constituent 4-simplices of 
the 4-block. Each coefficient iyA''®®- where X = 1, 2, 3, 
denotes the number of AX, BX triangular hinge pairs in 
a Cauchy surface and would equal the number of tetra¬ 
hedral edges that generate the pairs. So TVA''®'^ i is given 
by twice the number of parent edges Ni, that is. 


Ni 271 - n6»(2) „ 271 - 20(1) - 40(3) 

- , _ = 3- 




3/3 


271 - 40 ( 4 ) 




(141) 


VArea 1 = 2Ni. (138) 

As there are three Ui-edges per parent face, then iVA"'®a 2 
is 

VArea 2 = 3 ^^ = 6 V 3 , (139) 

where N 2 is the number of parent faces, and the second 
equality follows from (11321) . Finally, iVA"'®a 3 is deter¬ 
mined by the fact that each parent tetrahedron has six 
Pi-edges, none of which is shared with any other parent 
tetrahedra. Therefore, 

VArea 3 = 0 JV 3 . (140) 


where relation (11161) has been used to simplify the ex¬ 
pression. From (jl28|) . we can parametrise v in terms of 
cos 0 (^) through the expression 

.2 12 (a^ — I) [(4 q; 2 — l) cos 0 ( 2 ) — (2a‘^ — l)] 

q; 4 (l -p cos 0 ( 2 )) 

(142) 

As with the parent model, the range of 0(^) is bounded 
from above by the requirement that > 0 and from 
below by the requirement that the strut-lengths be time¬ 
like, that is, irrii)^ < 0 ; this leads to the range of 

^ < 0(2) < arccos • (143) 
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We have also varied the action with respect to Vi to 
obtain the evolution equation for the children models. 


The calculation, which we have not shown, is similar to 
that for the parent models. In the continuum time limit, 
we obtained 


0 = 


Ni 




{“-m) 


1/2 


Sna^vv^i) 


2a/q; 2 _ i 


3(4a2 - 1) - 




3- 




1/2 


16q; (27r — 




4 • 2 
a V 


4- i 




4 48(a2-i)" 2(a2-i) 4 


- 1 




1/2 


SVSvv^v 


4a2 (a2 I) + 2 (a2 - i) - l) - {a^ - |) - l) 




3(4a2 _ 1)_ 




8;5“ (;9= - I) + 4 (3= - 1) (/^T - l) - /3^ (/3^ -!)(?&- l) 'l 


2 _ i 

3 


4-(5ex-l)»4 


1/2 


1 _ (^!_|)_/;2 


48 (271-20(1)-40(3)) /i ^ 


,4„-,4 


a'^v'^ 


4- 




3 ( 4/32 - 1 ) - 


- -1 


I3^v^ 

4(/3=-5) 


_ 1 ^2j,2 


1/2 


3l3°vv‘^v 


'2 _ i 


13‘^v^ 


3(4/3^ - 1) - 

8/3(271-40(4)) / ^4 


3- 




4(/3^-3) 


1/2 




-nti — 


/34/>2 /3 


2^2 


4-4^ V'^ 48 (/32 -i)^ 2(/32-i) 

12-\/3 (a^ — i) a^ui) 24-\/3 (/3^ — i)/3^i;ii 


- 1 


a2-i 


3 (4a2 _ 1 ) _ 




2 _ 1 


3(4^2_i)_ 


, 3V3 .. 

+ —VV 


A 


/ 


a2 - ^ 
^ 3 


3^-1 


\ 


- i 1 - 


F-i), 




4^3 


+ 2. 


a2-i 


1 , F + 3h2)-12(^a2_ l+2^/32_i 


(144) 


If this evolution equation is to be consistent with the Hamiltonian constraint (I141D and with relation (I142L 
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then we require the Hamiltonian constraint to be its first 
integral. In Appendix |51 we have proven that this is pos¬ 
sible if a = /3 = 1, that is, if the children tetrahedra are 
all equilateral, and we have found indications that this 
requirement fails for any other values of a and /3. 

By taking a = /3 = 1, we can further simplify our 
model. The dihedral angles now become 

cos0^^^ = cos0^^^ = cos 9^^^ = cosO^'^^ 

_ 1 + (145) 

which is identical to its counterpart (H71) for the parent 
models. Since all dihedral angles are identical, we shall 
henceforth drop the superscript. Then, the Hamiltonian 
constraint (11411) simplifies to 


V 


2 


1 


Ni 

(27r — nO) + 6 (27r — 50) , 

As J 


(146) 

and relation (I142D for v simplifies to 


V 


2 


1 — 2 tan^ 



(147) 


which is identical to its parent model counterpart as given 
by (l56l) . We can now simply use (11461) and (11471) rather 
than the evolution equation to determine the behaviour 
of our models. 

Finally, we present the volume of the subdivided Regge 
universe, which is again simply the sum of the volumes 
of all constituent tetrahedra in a Cauchy surface. This 
volume is 


Un, ik) = ^ As (l + 2 ^ (148) 

which reduces to 


UN^[U) = V2N3v{Uf, (149) 


when a = fi = 1. 


E. Initial value equation for the children 
models 

Like the parent models, the children models also admit 
a moment of time symmetry at the moment of minimum 
expansion, when u = 0. At this moment, it follows from 
(11411) that 


vl = 


Ai 


-I [^ + 2 


3 


+ 3/3 ( 2 ^- 40 ^^)) 


27r - - 40^^^' 


where 


cos 0Q^^ 

cos 0Q^^ 
cos 9q^^ 

cos 9lf^ 


1 

V'3(4a2-1)’ 
2a^ - 1 
4a2 _ 1 ’ 

1 

V 3 ( 4 / 32 -l)’ 
2/32 - 1 
4/32 - 1 ■ 


The three deficit angles given, respectively, 

by (11311) . (11331) . and (11341) . now become 


= 271 — n9g^^ 

5 ( 2 ) = 27r - 20^^^ - 40^^^ (151) 

4 ^) = 27r-40^'‘\ 


We shall now demonstrate that (11501) is also consistent 
with the initial value equation. Unlike the parent mod¬ 
els however, vertices in children Cauchy surfaces are no 
longer identical, so the notion of a ‘volume per vertex’ be¬ 
comes harder to define; thus we cannot apply Wheeler’s 
definition of as given by (1571) . to the initial value 
equation (l84l) . Instead, we shall integrate (l84l) over the 
entire Cauchy surface Eq. The left-hand side becomes [l| 

f ^^'>RdPx = 2^ (.iSi, 

2 e {hinges} 


where the integration measure is unity because the Regge 
tetrahedra are flat, the summation is over all edges in Eq 
because these are the hinges of a 3-dimensional skeleton, 
5i is the 3-dimensional deficit angle of an edge, and is 
the edge’s length. The initial value equation (IMl) there¬ 
fore becomes 


iiCi = 8 tt p cPx. (152) 

i e {hinges} 

To the best of our knowledge, this form of the Regge 
initial value equation is novel and applies generally to any 
Regge Cauchy surface at a moment of time symmetry. 
For our particular model where the matter content is 
simply the cosmological constant A, once again we have 
that Pa = A/Stt, and hence, the integration on the right- 
hand side simplifies to AUN^{ti). Making use of (11481) . 
we can therefore express the initial value equation for the 
children models as 

2G{edges} V \ / 

(153) 

Recall that a child Cauchy surface has only three dis¬ 
tinct edge-lengths, which for surface Eq are Uq, Vq, and 
Po, and that all edges with the same length are iden¬ 
tical, implying that they would be associated with the 


( 150 ) 
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same deficit angle in the 3-dimensional skeleton of Eg. 
We shall denote the 3-dimensional deficit angle associ¬ 
ated with uo by 6q^\ with vo by Sq^\ and with po by 
Then the summation on the left-hand side of (11531) can 
be expressed as 

^ ^(1) + (^J^Edse j(2) 

iG{edges} 

where «o^ ]\^Edge po ^ ^nd denote the num¬ 

bers of edges with lengths uq, po, and vq, respectively, in 
Eg. Making use of the scaling relations (I89p and the fact 
that 


jyEdge Uo _ jyArea 1 
jyEdge Vo _ jyArea 2 
^Edge Po _ jyAiea 3 


we further simplify the summation to 

^ £A = 2 (W a + SNg + SNg p , 

2G{edges} 


where we have used relations (|138p to (11401) to substitute 
for \ iV^"'®®' 2, and iV^"®® 3. 

Substituting back into (I153D , we solve for ug to obtain 


Vn = 


2V3( 


A /32 


iV3^° 


+3^"^ +3/3^^^ 


(154) 


This expression is identical to (11501) provided the 3- 
dimensional deficit angles Sq'^ equal their 4-dimensional 
counterparts Sq \ as given by (I15ip . for i = 1, 2, 3. 

The deficit angles Sq^ and Sq^ have identical forms 
though. Recall that each edge in the hypersurface gener¬ 
ates a trapezoidal hinge, and each triangle meeting at the 
edge generates a 3-dimensional face meeting at the hinge. 
Thus from each of dehcit angles Sq^ , we can immediately 
deduce the form of its counterpart ; it follows that 

^ =2tt - n9Q^\ 

) = 27r - 20W _ 40(3)^ 

) = 2 ^ - 4^, 

where 9q^ is the 3-dimensional dihedral angle between 
the two triangles that generate the two 3-dimensional 
faces separated by dihedral angle 0 q . To complete our 
proof then, we need only demonstrate 9q^ = 9q^ for all 
i. 

Consider a typical tetrahedron in E3 with vertices 
A, B,C, D. Let the vertices’ co-ordinates be identical 


to the spatial co-ordinates of their counterparts in (lllll) . 
This tetrahedron has two distinct hinges, AB and AD. 
We hrst consider the dihedral angle at AB. Triangles 
ABC and ABD meet at this hinge and are separated 
by the dihedral angle 0g^\ The unit normal to ABC is 
(0,0,1), and the unit normal to ABD is 


1 



2 

0 






From their scalar product, we find that 


cos 0g^^ 


1 

V3(4a2-1)’ 


and thus we see that 9^'^ = ^g^^- 

Next, we consider the dihedral angle at hinge AD. Tri¬ 
angles ABD and ACD meet at this hinge and are sep- 

“( 2 ) 

arated by the dihedral angle 0g . The unit normal to 
ABD is the same as before, while the unit normal to 
ACD is 


1 

\/3 {4ul-vl) 





The scalar product gives 


COS 0 g^^ = 


2 a 2 _ 1 

4q;^ — 1 ’ 


~( 2 ) (2) 

and we can similarly conclude that 9q = 9q . 

Finally, we can readily obtain analogous expressions 
for 9q^^ and by swapping ug for po or equivalently 
a for /3 in all of the above expressions. From we 
immediately obtain 


cos 0g^^ 


1 

v'3(4/32-l)’ 


and from 9, 


(4) 
0 ’ 


we obtain 


cos 9q'^^ 


2^2 _ 1 
4/?2 - 1 ■ 


Therefore 0g^^ = 9q^^ and 9q'^^ = 9q'^\ and we can now 
conclude that the Regge equation (I14ip for the children 
models does indeed satisfy the initial value equation. 
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F. Discussion of the children models 

We now compare the behaviour of the children mod¬ 
els against their corresponding parent models as well as 
against the FLRW uiriverse. We shall only consider the 
models where the children tetrahedra are all equilateral, 
as these are the only models for which the evolution and 
constraint equations are consistent. 



U(t) 


Figure 13. The expansion rate of the universe’s volume 
against the volume itself for all parent models and their chil¬ 
dren models. 


Table V. The minimnm volume and the fractional difference 
from the FLRW minimum for each Regge model and for the 
FLRW model. 


Model 

Minimum volume 

Fractional 
difference from 
FLRW 

FLRW 

102.567937639753 

0 

5-tetrahedral 

parent 

171.741398309775 

0.67442 

16-tetrahedral 

parent 

135.70186007972 

0.32304 

600-tetrahedral 

parent 

105.692461545881 

0.03046 

60-tetrahedral 

child 

172.637934789289 

0.68316 

192-tetrahedral 

child 

179.191098180344 

0.74705 

7200-tetrahedral 

child 

1268.30953855058 

11.36556 

Figure US! and Figure [TH compare the graphs of dU/dt 


against U for each of our models, with U given by p49l) 
for the subdivided Regge models and by (150)) for the par¬ 
ent models. The graphs in Fig. [14] focus on the 600- 
tetrahedral parent and its 7200-tetrahedral child, with 
the bottom plot extended so as to reveal the graphs’ 
endpoints. Once again, the models reproduce the correct 
qualitative dynamics of the universe being approximated. 
However, at low volumes, subdividing the tetrahedra ac¬ 
tually made the accuracy of the Regge approximation 
worse. We can see this more concretely in Table fVl where 
we list the minimum volumes for each model and their 
fractional difference from the FLRW minimum. While 
increasing the number of tetrahedra in the parent mod¬ 
els brought the minimum volume closer to the FLRW 
value, increasing the number in the children models ac¬ 
tually brought it further away. In fact, the worst parent 
model was still more accurate than the best child model. 

We also see that all models again diverge from the 
FLRW model as the universe expanded; however increas¬ 
ing the number of tetrahedra reduces the rate of diver¬ 
gence, and in this sense, increasing the number of tetra¬ 
hedra improves the Regge approximation. We again be¬ 
lieve this lower rate of divergence is the result of more 
tetrahedra providing a higher resolution approximation 
that can ‘keep up’ longer with the FLRW model. We also 
note that each model terminates at large volumes when¬ 
ever the strut becomes null and that this endpoint gets 
increasingly delayed as the number of tetrahedra is in¬ 
creased. Thus these figures reveal that each child model 
provides an approximation that starts off worse than its 
parent but is later much better by virtue of its more 
robust resolution. Indeed, if one were to extrapolate 
all graphs past their end-points to very large volumes, 
the 7200-tetrahedral model would ultimately provide the 
best performance. 

As with the parent models, we can define a 3-sphere 
radius R{t) analogous to (IMll for children Cauchy sur¬ 
faces such that R{0) = a(0); that is, R{t) and a{t) match 
at the moment of minimum expansion; thus, we define 
R{t) to be 

i?(t) = —u(t), (155) 

'^min 

where Vmin is the minimum value of v and is given 
by (11411) when = arccos[(2a^ — l)/(4a^ — 1)] = 
arccos(l/3). Such a definition is possible because, as dis¬ 
cussed at the end of Section IIV A[ all dynamical length- 
scales in the system are related to each other by time- 
independent scalings such that there is really only one 
independent dynamical length-scale describing the entire 
model; therefore any dynamical length-scale, when re¬ 
scaled by an appropriate constant, will yield the same 
R{t). Note that when n = 5, then l{t) as given by (iTTl) 
and vit) as given by (11461) will be identical apart from 
an overall constant factor. When this happens, the two 
corresponding models will have the same R{t). This hap¬ 
pens for the 600-tetrahedral parent model and its 7200- 
tetrahedral child. 
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U{t) 

(a) 



U(t) 

(b) 


Figure 14. The expansion rate of the universe’s volume 
against the volume itself for the 600-tetrahedral parent and 
its 7200-tetrahedral child model; (a) focuses on the region 
around the origin while (b) shows both Regge graphs in their 
entirety. 

We now examine the behaviours of dR/dt versus R{t) 
and of the corresponding 3-sphere volumes. Figure 1151 
shows the relationship between dR/dt and R{t) for each 
of the Regge models and the relationship between da/dt 
and a{t) for the FLRW model. Figure fTOl shows the rela¬ 
tionship between the corresponding 3-sphere volumes and 
their rates of expansion. We note that in both figures, 
the graphs for the 600-tetrahedral model and its 7200- 
tetrahedral child coincide, as expected. The definition 
of R{t) has clearly removed any variability in the initial 
performance of the Regge models; now, the Regge models 
with more tetrahedra consistently outperform those with 
fewer. The rate of divergence from FLRW is again re- 



Figure 15. A combined graph of the radius’ expansion rate 
against the radius itself, using R{t) as the Regge radius. 

duced and the graphs’ end-points further delayed as the 
number of tetrahedra is increased. Thus in these graphs, 
increasing the number of tetrahedra clearly improves the 
Regge approximation. 

As mentioned at the start of Section llVDi we were un¬ 
able to obtain a consistent set of Regge equations when 
we locally varied the children models’ action. We sus¬ 
pect that we may not actually have complete freedom to 
specify the lengths of all three sets of struts, contrary 
to Brewin’s claim. Instead, two of the strut-lengths may 
depend on the third. The dependence may actually be 
determined by two of the three constraint equations ob¬ 
tained from varying the struts locally; they may perhaps 
act as evolution equations for the two non-independent 
sets of struts, one equation for one set; the third con¬ 
straint equation would still determine the evolution of 
the tetrahedral edges. We note though that so long as 
the strut-lengths remain interdependent, there will also 
be some constraint between the tetrahedral edge-length 
ratios ol and /3 similar to (11171) : this would follow from 
a similar reasoning to that which led to (11171) . After de¬ 
termining the lengths of the two non-independent sets of 
struts, one could then deduce from them the constraint 
between a and /3. 

Nevertheless, we may still run into the same prob¬ 
lem encountered when locally varying the parent model. 
There, we discovered that for the evolution equation (1751) 
to be a first integral of the Hamiltonian constraint (1591) . 
the momentum constraints dzil) must also be satisfied. 
However these latter constraints required the model to 
behave in an unphysical manner, and thus, the model 
broke down. We suspected the underlying cause to be 
the breaking of Copernican symmetries from introducing 
the diagonals. With the children models, we should also 
expect the diagonals to break Copernican symmetries, 
possibly rendering, for example, the mid-point vertices 
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(a) 



(b) 

Figure 16. A combined graph of the 3-sphere volume’s expan¬ 
sion rate against the volume itself, using R{t) as the Regge 
3-sphere radius; (a) focuses on the region around the origin 
while (b) shows all Regge graphs in their entirety. 


inequivalent to each other. Nevertheless, even if the local 
model were not viable, it would, through a chain-rule re¬ 
lationship analogous to (I76L produce an alternative but 
viable global model, and the properties of this new model 
would be worthy of further investigation. Indeed, it may 
possess some desirable advantages over our current child 
model, such as, for instance, having a Hamiltonian con¬ 
straint that is unconditionally a first integral of the evo¬ 
lution equation. However, we shall leave for now a more 
thorough examination into the viability and properties of 
such models, local and global, to future study. 
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APPENDIX A: AVERAGE RADIUS OF A 
CW CAUCHY SURFACE 

In this appendix, we shall present the computation of 
the numerical factors in (|28)) for the average radius of a 
CW Cauchy surface. To compute these factors, we must 
first determine the radius to any arbitrary point v in the 
Cauchy surface. In terms of the position of vertices A, 
B, C, D oi our representative tetrahedron, the position 
of V in the tetrahedron can be expressed as 

v = q:A-|-/IB-|-7C-|-(5D, (A1) 

where A, B, C, D are the position vectors of vertices A, 
B, C, D relative to the embedding 3-sphere’s centre, and 
where constants a, /3, 7 , S satisfy 0 < a, /?, 7 , 5 < 1 and 
a + P ^ + 5 = 1. Working in the E"^ co-ordinate system 
of (ITbl) . we can take vertices A, B, C, D to have co¬ 
ordinates given by vertices 1 to 4, respectively, in Table 
im Then the radius of v is given by 


J 


Rv{ti) = 


R{ti) 


1 — 2 cos ■ 


27r 


27r , 


1 — 2 cos — 1 -I- 2 cos — {aP+a'^+a5+P^+p5-\-^5) 


1/2 


(A2) 


where R{ti) is the radius to any of the vertices and is 
determined from the edge-lengths l{ti) using ( 1 ^ . 

To get the average radius, we therefore need to com¬ 
pute the multiple integral 

R{ti) = — / / / Rydadpd-f, (A3) 

J q;=0 J I3—0 J 7=0 


where the normalisation M is given by 

^1 nl — OL pl — OL — ^ -1 

J\f= / / dadp dj = -. (A4) 

Ja=0 Jp=0 J~/=0 b 

For each of the three models, the integral (IA3|) was eval¬ 
uated numerically to obtain the factors in (|28)) . 
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APPENDIX B: PROOF OF CONSTRAINT 
EQUATION BEING A FIRST INTEGRAL OF 
THE EVOLUTION EQUATION 

In this appendix, we shall first prove that the Hamil¬ 
tonian constraint (1571) is a first integral of the evolution 
equation (1551) for the parent models. We shall next pro¬ 
vide a partial proof that the Hamiltonian constraint (11411) 
is a first integral of the evolution equation (11441) for the 
children models. 

To facilitate comparison with the constraint equation, 
we shall first simplify the parent model’s evolution equa¬ 
tion. Note that the number of edges Ni in a Cauchy 
surface is related to the number of triangles N 2 by the 
relation Ni = 37 V 2 /n, where n is the number of trian¬ 
gles meeting at an edge; this relation follows because a 
triangle has three edges, but an edge joins n triangles 
together. By using this relation to substitute for N 2 and 
by substituting (1571) into the P factor in the last term of 


J 


I = 3V2 


Ni 1 



We next multiply through by I and replace the on the 
left-hand side by dsg. After further simplification, we 
arrive at 


1 {2tt — nO) 




n ll 

2 V 2 _ 1 


which is clearly identical to dHU apart from an irrelevant 
overall factor of (|Z^ — 1)/(1 — The numerator of 

I 


(l58)) . we can simplify (1581) to get 



Now we consider the constraint equation (1571) itself. To 
demonstrate that this is a first integral of (l58l) . we first 
differentiate it; this gives 


li = 3V2 


Ni 


N.A 


(1 - jp) 


1/2 


8 1 - U2 


(B2) 

The derivative 6 can be determined simply by differenti¬ 
ating 9 = arccos , yielding 


e = - 


ll 


2s/2 


(i-Jp) 


1/23-4/2’ 


and (IB2I) now becomes 


1 /■ 


1 - 



1/2 3 - i/2 


1 ;• (27r — nO) 

+ - i 

8 I-IP 


(B3) 


this factor is simply the square of the strut-length and is 
therefore constrained to be negative. This in turn implies 
that the denominator must be strictly positive. Hence, 
the factor is never singular nor zero. Therefore, we can 
conclude that the constraint equation is indeed a first 
integral of the evolution equation for the parent models. 

We now turn to the children models. Following the 
example of the parent model, we similarly differentiate 
constraint equation (jl41l) with respect to time and sim¬ 
plify the resulting expression to get 


0 = 




N. 


(a 1 


1/2 


a 2 -i 


3na^vv 


3(4a2 _ 1) _ 




3- 




1/2 


a (27r — 710 ^ 2 )^ /1 


4- 
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3 q;2 _ 
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18/3^?;?; 


3 ( 4/32 - 1 ) - 


4(/32-i) 


3- 


/3(27r-46l(4)) / /3 


4(/32-i) 


1/2 


1 

3 /32-1 


2/32i;2 


- 24 


(B4) 


If we use (I14ip to expand out the term at the end of 
(I144|) . we hope to find that this expanded form of (11441) 
will be identical to (IB4I) apart from some overall factor. 
If this is the case, then we shall have demonstrated that 
(11411) is a first integral of (I144p . 

Note that (|144p and (IB4I) can be separated 
into four distinct components corresponding to a 


{N 1 /N 3 ) (27r — term, a (27r — 20 ^) _ 46/(3)^ term, 

a (27r — 40*^^^) term, and everything else. To identify 
the factor, our approach will be to match each com¬ 
ponent independently of the rest. We begin with the 
{N 1 /N 3 ) { 2 tt — component which requires the fol¬ 

lowing equality to hold: 


16 



4 

■jVV + 

3 


1 


(a^ + 

2 - i) 


a^vv^v 

48 



factor X 


/I a'^ 


1 

3 


vv + — 


3 Q;2 _ 



The highest order terms on the left-hand side are the v'^ 
and vif^v terms, which are of order 172 higher than the 
right-hand side. Thus we make the ansatz that the factor 
is of the form 

factor = A + B 

for some constants A and B to be determined. By match¬ 
ing the numerical constants on the two sides, we find that 
A = —4. However if we now match the vv terms, we ob¬ 
tain the expression 

4 . 4 ^ _ _ 4 ^ 


This equality can only be true if a = ±1 or ± l/-\/3. 


However since a is the ratio between the lengths of two 
edges, we require a > 0 , and for the strut-length to be 
finite, we also require a > l/-\/3- Therefore the only 
acceptable solution is a = 1. After making this choice, 
we can continue matching the remaining terms to find 
that B = 3/2. Thus the factor is given by 

factor = —4 ^1 — ■ 

By matching the (27r — component in a similar 

manner, we also find that /3 = 1 and that the terms from 
the two equations also differ by the same factor. 

If we now make the substitution a = /3 = 1, we find 
that the evolution equation (11441) is indeed identical to 
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the time-derivative of the constraint equation (IB4I) mul¬ 
tiplied by the overall factor of —4[1 — ( 3 / 8 ) 1 ;^]. Therefore 


when a = /3 = 1, the constraint equation is a first integral 
of the evolution equation. 
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